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1.  INTRODUCTION: 


Most  prostate  cancers  can  be  managed  by  active  surveillance.  However,  the  high-risk  population, 
comprising  about  15%  of  prostate  cancer  cases,  has  a  significantly  higher  chance  of  local  recurrence  or 
distal  metastasis;  hence,  more  intensive  treatment  for  this  population  is  required.  Treatment  options,  such  as 
radiation,  surgery,  androgen-deprivation  therapy,  or  their  combinations,  are  often  ineffective  in  this  subset 
of  patients.  Disease  usually  relapses  after  therapy  and  there  are  no  effective  preventive  measures. 
Immunotherapy  holds  great  potential  in  treating  cancer  but  major  challenges  remain.  We  have  developed 
Tumor  Immunotherapy  by  Gene-circuit  Recruited  Immunomodulatory  Systems  (TIGRIS)  as  a  novel  cancer¬ 
detecting  gene  circuit  platform  that  enables  tumor-localized  immunomodulation.  The  outputs  of  the  circuit 
include  surface  T-cell  engagers  (STEs),  which  can  trigger  the  antigen-independent  T-cell  killing  of  tumor 
cells.  TIGRIS  is  a  programmable  and  highly  specific  platform  for  triggering  cancer  immunotherapies  that 
can  be  easily  adapted  to  various  types  of  cancer,  such  as  prostate  cancer.  Cancer-detecting  gene  circuits  will 
trigger  T  cells  to  kill  tumor  cells  but  leave  normal  cells  unharmed.  Immune  memory  may  further  develop 
after  TIGRIS  administration  and  may  prevent  future  cancer  relapse.  If  successful,  TIGRIS  would  obviate  the 
need  for  traditional  therapies,  such  as  chemotherapies,  that  have  high  systemic  toxicity,  and  would  improve 
the  survival  and  life  quality  of  prostate  cancer  patients. 


2.  KEYWORDS:  cancer-detecting  gene  circuit;  tumor  immunotherapy;  prostate  cancer;  metastasis 

3.  ACCOMPLISHMENTS: 

What  were  the  major  goals  of  the  project? 


Training-Specific  Tasks: 


Major  Task  1:  Training  and  educational  development 
in  synthetic  biology 

Months 

Percentage  of 
Completion 

Subtask  1 :  Attend  relevant  seminars  at  Synthetic 

Biology  Center  and  Broad  Institute  at  MIT 

1-24 

50% 

Subtask  2:  Present  research  at  the  weekly  Lu  lab 
focused-group  meetings 

1-24 

50% 

Subtask  3 :  Attend  a  national  scientific  meeting  in 
synthetic  biology 

1-24 

50% 

Milestone(s)  Achieved:  Presentation  of  project  data  at  a 
national  meeting  or  publishing  of  project  result  in 
relevant  journals 

24 

Reached 

Major  Task  2:  Training  and  educational 
development  in  prostate  cancer  research 

Months 

Percentage  of 
Completion 

Subtask  1 :  Attend  relevant  seminars  at  Koch  institute  for 
Integrative  Cancer  Research  at  MIT 

1-24 

50% 

4 


Subtask  2:  Attend  a  national  scientific  meeting  in 
prostate  cancer  research 

1-24 

0% 

Milestone(s)  Achieved:  Presentation  of  project  data  at  a 
national  cancer  related  meeting  or  publishing  of  project 
result  in  relevant  journals 

24 

Not  yet  reached 

Research-Specific  Tasks: 


Specific  Aim  1 :  Design  and  build  RNA-based  AND  gate 
for  Tumor  Immunotherapy  by  Gene-circuit  Recruited 
Immunomodulatory  Systems  (TIGRIS) 

Timeline 

Percentage  of 
Completion 

Major  Task  1:  Create  and  optimize  the  genetic 
components  for  building  a  high  performance  single¬ 
output  RNA-based  AND  gate  in  mammalian  cells 

Months 

Percentage  of 
Completion 

Subtask  1 :  Design  and  create  prototypes  of  the  RNA-based 
AND  gate  in  mammalian  cells. 

1-6 

85% 

Subtask  2:  Optimize  each  component  of  the  RNA-based 
AND  gate  to  achieve  high  ON: OFF  performance. 

3-9 

100% 

Milestone(s)  Achieved:  Develop  a  modular  single-output 
RNA-based  AND  gate  that  achieve  high  ON: OFF 
performance  in  mammalian  cells. 

9 

Not  yet  reached 

Major  Task  2:  Extend  the  single-output  RNA-based 

AND  gate  to  trigger  multiple  outputs. 

Months 

Percentage  of 
Completion 

Subtask  1 :  Implement  an  additional  layer  of  transcription 
cascade  to  enable  RNA-based  AND  gate  to  express 
multiple  outputs. 

9-15 

0% 

Subtask  2:  Optimize  the  AND-gate  topology  to  enhance 
the  ON: OFF  performance. 

9-15 

0% 

Milestone(s)  Achieved:  Develop  a  multi-output  AND  gate 
with  high  ON: OFF  performance. 

15 

Not  yet  reached 

Specific  Aim  2:  Validate  whether  the  multi-output 

AND  gate  can  be  adapted  to  distinguish  cancer  cells 
from  normal  cells. 

Major  Task  3:  Characterize  the  output  level  of  AND 
gate  in  tumor  cells  and  normal  cells 

Months 

Percentage  of 
Completion 

Subtask  1 :  Identify  ideal  promoters  that  can  sense  prostate- 
cancer- specific  transcriptional  activity 

9-15 

30% 

Subtask  2:  Build  prostate-cancer-specific  AND  gate  that 
can  distinguish  prostate  cancer  cells  from  normal  cells. 

15-24 

0% 

Milestone(s)  Achieved:  Demonstrate  and  optimize  a  multi¬ 
output  AND  gate  that  can  distinguish  prostate  cancer  cells 
from  normal  cells. 

24 

Not  yet  reached 

What  was  accomplished  under  these  goals? 

Training-specific  tasks: 

I  (Ming-Ru  Wu)  have  been  actively  attending  seminars  at  the  Synthetic  Biology  Center,  Dana-Farber 
Cancer  Institute,  Broad  Institute,  and  Koch  Institute  for  Integrative  Cancer  Research  at  MIT.  Seminars 
covered  topics  major  topics  such  as  cancer  biology,  cancer  genomics,  tumor  immunology,  and  immune 
engineering.  I  gained  a  lot  of  expertise  and  am  up  to  date  on  the  current  advances  in  cancer  therapy.  I  have 
also  attended  an  international  conference  on  synthetic  biology,  the  4th  International  Mammalian  Synthetic 
Biology  Workshop  at  Boston  this  year  (2017/5/6-2017/5/7).  I  also  attended  an  international  conference  on 
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bioinformatics,  Bioconductor  BioC2017  conference  at  Boston  (2017/7/27-2017/7/28).  I  have  presented  this 
work  multiple  times  in  lab  meeting  or  local  biological  meetups  during  this  training  period.  I  have  also 
published  our  work  as  co-first  author  in  the  highly  ranked  journal  Cell  (see  Appendix  for  more  detail).  I  will 
be  giving  a  talk  on  this  project  at  the  next  Annual  Symposium  of  Synthetic  Biology  Center  at  MIT  (January 
2018). 

In  addition  to  the  above  conference  and  lab  meeting  training  activities,  I  have  also  obtained  a  tremendous 
understanding  of  synthetic  biology,  both  theoretical  and  practical.  Furthermore,  I  have  gained  extensive 
hands-on  experience  in  designing,  building,  and  validating  cancer-detecting  gene  circuits.  I  have  also  gained 
hands-on-experience  on  establishing  advanced  tumor  models,  such  as  establishing  an  advanced  xenograft 
model  in  NSG  mice. 

Research-specific  tasks: 

One  of  the  major  goals  during  this  reporting  period  has  been  to  validate  the  hypothesis  that  cancer¬ 
detecting  gene  circuits  that  trigger  immunotherapy  can  be  built.  Most  gene  circuits  utilize  protein  molecules 
for  computation.  However,  those  protein  molecules  may  be  immunogenic  for  long-term  clinical 
applications.  Therefore,  we  focused  on  building  circuits  that  only  use  RNA  molecules  for  computation.  If 
successful,  these  circuits  will  be  a  major  advancement  in  synthetic  biology  and  in  the  future  treatment  of 
cancer  patients.  We  first  designed  a  novel  architecture  to  develop  RNA-only  AND  gate  computation  (Figure 
1).  To  ensure  that  these  circuits  perform  as  intended,  we  systematically  optimized  each  genetic  component 
that  serves  as  a  basic  building  block  for  the  circuits.  For  example,  we  have  optimized  the  microRNA 
(miRNA)  binding  site  sequences,  miRNA  sponge  architectures,  and  synthetic  intron  architectures  for 
reaching  a  high  ON: OFF  ratio  (Figure  2-4).  We  then  utilized  these  components  to  build  the  AND  gate 
computation  unit.  We  further  enabled  the  circuit  to  produce  multiple  outputs  by  utilizing  the  synthetic 
transcriptional  factor  GAD  (a  fusion  protein  consisting  of  the  GAL4  DNA  binding  domain  and  the  VP  16 
transcription  activating  domain)  as  the  output  of  the  first  layer  of  the  circuit.  GAD  can  then  further  trigger 
the  second  layer  of  the  circuit  and  thereby  enable  the  production  of  multiple  outputs  (Figure  5A).  We  have 
validated  that  this  circuit  can  distinguish  tumor  cells  from  normal  cells  (Figure  5B),  and  indeed,  it  enhances 
anti-tumor  specificity  (Figure  5C).  This  optimization  process  was  done  with  an  ovarian  cancer  model; 
nevertheless,  we  can  adapt  the  circuit  to  target  other  types  of  cancer  cells.  We  have  demonstrated  this 
adaptability  by  switching  the  sensors  to  breast  cancer-specific  sensors  (Figure  5D).  This  observation 
convinced  us  that  we  can  readily  adapt  this  circuit  to  treat  prostate  cancer,  as  long  as  we  can  find  prostate 
cancer  sensors.  We  have  also  found  sensors  for  detecting  prostate  cancers  (will  be  discussed  in  detail  in  the 
following  paragraphs).  Furthermore,  we  have  validated  that  TIGRIS  has  the  potential  to  trigger  robust 
immunotherapy  in  vitro  and  in  vivo  (Figure  6-7).  Our  circuits  trigger  T  cells  to  kill  tumor  cells  and  secrete 
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IFN-y  specifically  (Figure  6A-C).  Our  circuits  can  also  trigger  tumor-specific  production  of 
immunomodulators,  such  as  STE,  CCL21,  IL12,  and  anti-PDl  antibody  in  vitro  (Figure  6D-G).  These 
immunomodulators  are  also  effective  for  treating  prostate  cancer.  Furthermore,  these  circuits  trigger  robust 
therapeutic  efficacy  in  vivo  even  when  only  15%  of  the  tumor  cells  bear  the  circuit  (Figure  7). 

We  have  done  further  work  on  adapting  this  circuit  to  treat  prostate  cancer.  By  collaboration  with  other 
postdocs  in  the  lab,  we  have  built  a  library  of  synthetic  tumor  sensors  (  >  6000  different  synthetic 
promoters)  for  the  TIGRIS  system.  This  library  can  be  used  to  detect  multiple  different  cancer  types. 
Currently,  we  are  looking  for  new  sensors  to  develop  this  system  for  the  treatment  of  prostate  cancer  (Table 
1).  When  14  breast-cancer  specific  promoters  were  tested  on  prostate  cancer  cell  lines  PC3  and  DU145,  two 
of  them  exhibited  activity  that  was  more  than  nine-fold  higher  than  background.  We  are  using  these 
promoters  to  develop  immunomodulatory  gene  circuits  specific  to  prostate  cancer. 

In  summary,  we  have  optimized  the  key  circuit  building  blocks  for  creating  highly  effective  RNA-based 
synthetic  gene  circuits  for  tumor-specific  immunomodulation.  Thus,  we  have  put  in  place  the  basic  structure 
of  a  system  that  is  theoretically  applicable  to  many  types  of  cancer.  We  have  also  optimized  circuit 
architecture  to  be  highly  modular,  thus  readily  adaptable  to  other  cancer  types,  such  as  prostrate  cancer,  and 
future  clinical  translation.  We  have  also  built  a  sensor  library  consisting  of  synthetic  sensors  that  perform 
better  than  native  sensors  for  targeting  cancers  (Figure  8)  and  that  serve  as  tumor  sensors  for  the  TIGRIS 
system.  The  next  step  forward  is  to  validate  the  specificity  of  TIGRIS  for  prostate  cancer.  Once  the  sensors 
have  been  combined  with  an  optimized  circuit  architecture,  we  will  determine  whether  the  system  remains 
highly  specific  to  prostate  cancer  cells  but  does  not  affect  other  major  tissue  types  in  vitro. 

What  opportunities  for  training  and  professional  development  has  the  project  provided? 

The  conferences  and  seminars  I  attended  have  familiarized  me  with  the  latest  advances  in  the  general 
field  of  cancer  biology.  The  project  has  provided  me  with  extensive  hands-on  experience  on  designing 
cancer-detecting  gene  circuits,  optimizing  circuit  inputs  and  outputs,  and  cancer  immunotherapy.  It  has  also 
provided  me  with  hands-on  experience  in  working  with  prostate  cancer  models. 

How  were  the  results  disseminated  to  communities  of  interest? 

Our  work  drew  significant  attention  from  the  cancer  therapy  field.  Our  work  was  highlighted  by  MIT 
news  by  journals  such  as  Cancer  Discovery,  Nature  Reviews  Cancer,  and  Nature  Reviews  Drug  Discovery. 
For  your  convenience  and  information,  please  see  the  following  links  for  details. 

Featured  in  MIT  news: 

http://news.mit.edu/2017/gene-circuit-switches-inside-cancer-cells-triggers-immune-attack-1019 
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Highlighted  in  Research  Watch  by  Cancer  Discovery  Journal: 

http://cancerdiscoverv.aacrioumals.org/content/early/20 1 7/10/27/2 1 59-8290. CD-RW20 17-20 1 

Highlighted  in  Research  Highlight  by  Nature  Reviews  Drug  Discovery: 
https://www.nature.com/articles/nrd.20 1 7.242 

Highlighted  in  Research  Highlight  by  Nature  Reviews  Cancer: 
https://www.nature.com/articles/nrc.20 17.11 2?WT.mc  id=TWT  NatureRevCancer 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

We  will  utilize  the  identified  sensors  to  build  TIGRIS  system  for  prostate  cancer.  First,  we  will  validate 
that  the  system  is  specific  to  prostate  cancer  cells  but  that  it  does  not  affect  normal  prostate  cells  or  vital 
normal  tissues.  Second,  we  will  further  optimize  the  circuit  to  minimize  the  potential  off-target  effect  on  a 
panel  of  normal  cells.  Lastly,  we  will  identify  the  optimal  circuit  architecture  for  reliably  triggering  multiple 
output  production.  This  will  pave  the  way  to  achieving  successful  translation  for  building  therapeutic 
relevant  circuit  for  treating  prostate  cancers. 
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Figure  1.  The  architecture  of  basic  RNA-only  single-output  AND  gate. 

(A)  Module  1  of  the  AND  gate  is  built  as  an  auto-inhibitory  loop  such  that  it  represses  its  own  output.  Module  2  is 
designed  to  inhibit  the  auto-inhibition  of  Module  1.  Cancer-specific  promoters  PI  and  P2,  are  used  to  regulate  Module 
1  and  Module  2,  respectively.  The  output  from  Module  1  is  expressed  at  a  high  level  only  when  both  PI  and  P2  are 
active,  which  enhances  the  tumor  specificity  of  the  circuit  than  targeting  tumor  with  a  single  promoter. 

(B)  All  4  possible  input  states  and  their  respective  output  states  are  shown  for  the  RNA-only  single-output  AND  gate. 
Input  states  are  defined  within  the  square  brackets  by  whether  Module  1  and  Module  2  are  active,  where  0  means 
inactive  and  1  means  active.  In  State  [1,0],  PI  is  active  and  the  mKate2  transcript  is  expressed  (mK-Exl  and  mK-Ex2 
denote  mKate2  exonl  and  mKate2  exon2,  respectively).  However,  this  transcript  encodes  an  intronic  synthetic 
miRNA  (miRI)  within  the  mKate2  gene  that  inhibits  the  mKate2  transcript  by  targeting  perfect  match  miRI  binding 
sites  (BS(Pe))  in  the  3’  untranslated  region.  Thus,  the  mKate2  levels  are  minimal.  In  State  [0,1],  PI  is  inactive  so  the 
output  protein  mKate2  is  not  expressed.  In  State  [0,0],  neither  PI  nor  P2  are  active  and  thus  the  output  protein  mKate2 
is  not  expressed.  In  State  [1,1],  the  P2  promoter  expresses  a  sponge  for  miRI  that  is  based  on  bulged  miRI  binding 
sites  (BS(B)).  This  enables  sequestration  of  miRI  away  from  inhibiting  the  mKate2  transcript  expressed  by  PI,  thus 
allowing  for  mKate2  expression. 
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2>fdd-acti  vation>  1 J5 

Bv4  5'  GGGGTTGAAG^S^TAATTAAA  3' 

Bvl  5'  CCGCTTGAAG§TCTTTAATTAAA  3' 
Bv14  5'  GGGCTTGAAGAjefgg^ATTAAA  3' 

Bvl  1  5'  CCGCTTGAAA^ggTAATTAAA.  3' 

Pe  5'  CCGCTTGEWVGTCT^  3' 

Bv6  5'  CCGCTTGAAG^TTTAATTAAA  3' 

1  .5>  fold-activation 

Bvl  2  5'  CCGCTTGAAA^gJj^AATTAAA  3' 

Bvl  5  5'  CCGCTTGAAG^KcTTTAATTAAA  3' 

Bv9  5'  CCGCTTGAAWfleiTTAATTAAA  3' 
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Figure  2.  Optimization  of  the  miRNA  binding-site  sequences  in  the  sponge  in  Module  2. 

(A)  The  activity  of  the  promoters  SSXIp,  H2Alp,  and  CMVp  in  HEK-293T  cells  that  were  used  for  circuit 
optimization. 

(B) -(C)  Various  miRNA  binding  sequences  (miRI-BSs)  used  in  the  sponge  affect  sponging  activity,  assayed  by 
mKate2  fold-activation  by  each  sponge  (state  [1,1]  versus  state  [1,0]  in  Figure  1).  Nucleotides  highlighted  in  dark  blue 
are  mismatched  nucleotides.  Nucleotides  highlighted  in  red  are  insertions.  Light  blue  capital-delta  (A)  indicates 
deleted  nucleotides.  SSXIp  was  used  to  drive  Module  1  expression  and  H2Alp  was  used  to  drive  Module  2 
expression.  Module  1  encoded  a  transcript  with  miRl  as  an  intron  within  the  mKate2  gene,  where  miRl  targeted  3 
perfect  match  miRl-BS  encoded  downstream  of  mKate2.  Module  2  encoded  the  miRl  sponges,  described  above, 
downstream  of  an  ECFP  gene. 

(D)  MFI  of  mKate2  and  ECFP  in  binding  site  optimization  experiments.  Error  bars  represent  S.E.M.,  N  =  3  biological 
replicates.  MFI  =  Mean  Fluorescence  Intensity. 
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Figure  3.  Optimization  of  the  sponge  architecture  in  Module  2. 

(A)  The  number  of  miRNA  binding  sites,  the  sponge  expression  level,  and  the  architecture  of  the  sponge  affect 
sponging  activity.  Fold  activation  was  assayed  by  mKate2  fold-activation  by  each  sponge  (state  [1,1]  versus  state 
[1,0]  in  Figure  1).  SSXlp  was  used  to  drive  Module  1  expression  and  H2Alp  or  CMVp  was  used  to  drive  Module  2 
expression,  as  indicated.  Module  1  encoded  a  transcript  with  miRl  as  an  intron  within  the  mKate2  gene,  where  miRl 
targeted  3  perfect  match  miRl-BS  encoded  downstream  of  mKate2.  Module  2  encoded  the  miRl  sponges.  5B  and 
10B  refers  to  5  or  10  bulged  miRl  binding  sites  in  the  sponge  based  on  the  Bv2  design  from  Figure  SIB.  EC1  and 
EC2  correspond  to  exon  1  and  exon  2  of  ECFP  and  constructs  in  which  10B  was  encoded  within  the  ECFP  transcript 
with  different  architectures  (Cons  and  ZRANB1).  12TuD  (Tough  Decoy)  architecture  refers  to  24  miRl-BSs  encoded 
within  12  loops  of  an  RNA  hairpin  (2  miRl-BSs  per  hairpin  loop).  5B/100  bp  refers  to  5  bulged  miRl-BSs  separated 
by  lOObp  linkers  rather  than  5bp  in  all  other  sponges.  Error  bars  represent  S.E.M.,  N  =  3  biological  replicates. 

(B)  MFI  of  mKate2  and  ECFP  in  sponge  architecture  optimization  experiments. 
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Figure  4.  pri-miRNA  backbones  used  to  produce  miRI  from  Module  1  affect  circuit  performance. 

(A)  Module  1  variants  Mvl-Mv3  expressed  a  transcript  in  which  the  different  miRI  pri-miRNA  architectures  shown 
were  encoded  as  an  intron  within  the  mKate2  gene,  where  miRI  then  targeted  3  perfect  match  miRl-BS  encoded 
downstream  of  mKate2.  The  colored  circles  denote  the  location  of  the  insertion  sequences  for  creating  Mv2  and  Mv3. 
(B-E)  Different  pri-miRNA  backbones  modulate  the  production  of  miRNA  and  the  ON:OFF  ratio  of  the  circuit.  (B-C) 
We  tested  the  three  different  Module  1  variants  together  with  Module  2  variations  composed  of  a  control  sponge  with 
no  miRI  binding  sites  (‘Cont.’)  and  a  sponge  with  5  perfect  match  miRI  binding  sites  downstream  to  an  ECFP  gene 
(‘5Pe’).  Module  1  was  expressed  by  SSXlp  and  Module  2  was  expressed  by  H2Alp.  Mvl-Mv3  exhibited  different 
miRI  production  efficiencies,  as  measured  by  ECFP  signal  reduction  with  the  5Pe  Module  2  compared  to  ECFP  signal 
in  the  control  Module  2;  Mv3  had  the  most  efficient  miRI  production. 

(D-E)  mKate2  expression  levels  with  control  sponge  or  the  optimized  sponge  (‘Sponge  7’,  Sponge  7  from  Figure  3). 
Mv3  exhibited  a  low  mKate2  level  in  the  OFF  state  and  the  highest  fold-activation  by  the  sponge  in  the  ON  state. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates. 
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Figure  5.  Multi-output  AND  gate  architecture  and  performance. 

(A)  The  incorporation  of  a  synthetic  transcription  factor  (GAD)  as  the  output  of  Module  1  means  that  it  is  only 
expressed  at  high  levels  when  both  PI  and  P2  are  high  (state  [1,1]).  GAD  enables  facile  tuning  of  downstream  gene 
expression  (mKate2  •  •  •  Output  N). 

(B)  The  G8-Pe  AND  gate  architecture  triggered  high  mKate2  expression  in  human  ovarian  cancer  cells  (OV8)  but  not 
in  normal  cells  in  state  [1,1]. 

(C)  Our  circuit  design  prevents  potential  off-target  effects  in  primary  human  T  cells.  The  S(E2Fl)p  promoter  was 
active  in  T  cells  but  S(cMYC)p  was  inactive  in  T  cells.  The  G8-Pe  circuit  triggered  minimal  mKate2  output  in  T  cells, 
thus  minimizing  the  potential  off-target  effects  that  would  have  been  observed  if  the  S(E2Fl)p  was  used  to  directly 
drive  mKate2  expression. 

(D)  Our  circuit  design  can  be  readily  adapted  to  distinguish  a  different  tumor  type  from  its  normal  counterparts.  By 
using  two  different  promoters,  the  G8-Pe  circuit  triggered  a  strong  output  in  the  breast  cancer  cell  line  MDA-MB-453 
but  minimal  output  in  the  non-tumorigenic  breast  epithelial  cell  line  MCF-10A. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates  for  all  experiments  except  N  =  6  biological  replicates  for 
OVCAR8  group  in  5B. 
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Figure  6.  Synthetic  circuits  triggered  human  ovarian  cancer-specific  expression  of  STE,  T-cell  mediated  lysis  and  IFN-y  secretion 
on  tumor  cells.  aHDF  (adult  human  dermal  fibroblast)  and  HOV  (human  ovarian  epithelium)  are  normal  primary  cells. 

(A)-(C)  The  circuit  in  Figure  5 A  was  engineered  to  display  STE  on  cell  surfaces  as  the  output.  STE  was  encoded  in  Module  3. 
State  4  [1,1]  cont.’  indicates  cells  containing  Module  1  and  Module  2  and  a  negative  control  Module  3 con  that  expresses  the  non¬ 
specific  transcription  factor  rtTA3  as  the  output.  State  [1,1]  indicates  cells  containing  Module  1,  Module  2,  and  Module  3.  State 
[1,0]  indicates  cells  containing  Module  1,  a  negative  control  Module  2con  in  which  S(E2Fl)p  expresses  ECFP  without  miRl- 
BS(B),  and  Module  3.  Circuit-transduced  tumor  cells  but  not  normal  cells  expressed  high  levels  of  STE  in  state  [1,1],  which 
triggered  robust  T-cell  killing  and  IFN-y  secretion.  Student’s  t-test  was  performed  to  compare  the  [1,1]  states  between  normal  cells 
and  OV8  cells,  and  the  [1,1]  state  versus  the  [1,0]  state  in  OV8  cells. 

(D)-(G)  Synthetic  circuits  triggered  human  ovarian  cancer-specific  expression  of  combinatorial  immunomodulators.  The  circuit  in 
Figure  5 A  was  extended  to  express  multiple  immunomodulatory  outputs  (surface-displayed  STE  and  secreted  CCF21,  IF12,  and 
anti-PDI  Ab)  specifically  in  OV8  cells  but  not  normal  primary  cells.  State  4 [1,1]  cont.’  indicates  cells  containing  Module  1  and 
Module  2  and  a  negative  control  Module  3con  that  expresses  the  non-specific  transcription  factor  rtTA3  as  the  output.  State  [1,1] 
indicates  cells  containing  Module  1,  Module  2,  and  Module  3.  Each  immunomodulatory  gene  was  expressed  from  its  own  G8p 
promoter  encoded  on  a  lentivirus,  except  for  the  anti-PDI  Ab,  which  was  split  into  two  lentiviruses  encoding  light  chain  (EC)  and 
heavy  chain  (HC)  respectively.  Cell  lines  were  co-infected  with  these  lentiviral  constructs.  Student’s  t-test  was  performed  to 
compare  output  levels  in  the  [1,1]  states  between  normal  cells  and  OV8  cells. 

Error  bars  represent  S.E.M.,  N  =  3  biological  replicates  (*  p  <  0.05;  **  p<0.005). 
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Figure  7.  Combination  immunotherapies  triggered  by  cancer-specific  circuits  significantly  reduced  tumor  burden  and  increased 
survival  in  an  intraperitoneally  disseminated  human  ovarian  cancer  model  in  mice. 

(A)  OV8  cells  were  transduced  with  constructs  encoding  a  control  output  rtTA3  (Control),  STE-only  (STE),  or  combination 
immunomodulators  (SCIP).  All  groups  were  implanted  into  NSG  mice  and  injected  periodically  with  human  T  cells.  Student’s  t- 
test  was  performed  to  compare  tumor  burden  between  groups  at  day  37.  Error  bars  represent  S.E.M.,  N  =  5  biological  replicates. 

(B)  Kaplan-Meier  survival  curves  of  various  groups  (N  =  5  biological  replicates).  A  log-rank  (Mantel-Cox)  test  was  performed  to 
compare  survival  between  groups. 

(C)  Robust  therapeutic  efficacy  was  achieved  even  when  only  15%  or  30%  of  tumor  cells  were  transduced  with  the  SCIP  circuit 
(the  remaining  implanted  cancer  cells  were  wild- type  OV8).  Therapeutic  efficacy  was  dependent  on  the  presence  of  T  cells. 
Student’s  t-test  was  performed  to  compare  tumor  burden  between  groups  at  day  41.  Error  bars  represent  S.E.M.,  N  =  4  biological 
replicates. 

(D)  Kaplan-Meier  survival  curves  of  various  groups  with  different  percentages  of  the  overall  tumor  cells  transduced  with  the 
SCIP  circuit  (N  =  4  biological  replicates).  A  log-rank  (Mantel-Cox)  test  was  performed  to  compare  survival  between  groups. 

(E)  Lentiviral  delivery  of  the  SCIP-expressing  circuit  significantly  reduced  ovarian  cancer  burden.  NSG  mice  were  injected  with 
OVCAR8  cells  into  the  peritoneal  space  at  day  0.  Lentiviruses  encoding  circuit  Module  1  and  Module  2  and  a  control  output 
rtTA3  (Control),  or  combination  immunomodulators  (SCIP)  were  then  injected  i.p.  on  day  7.  These  mice  were  also  injected  i.p. 
periodically  (day  9,  16,  and  23)  with  human  T  cells.  Student’s  t-test  was  performed  to  compare  tumor  burden  between  groups  at 
day  39.  Error  bars  represent  S.E.M.,  N  =  5  biological  replicates. 

(F)  Kaplan-Meier  survival  curves  of  various  groups  (N  =  5  biological  replicates).  A  log-rank  (Mantel-Cox)  test  was  performed  to 
compare  survival  between  groups.  (*  p  <  0.05;  **  p<0.005). 
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Figure  8  -  Specificity  comparison  of  natural  and  synthetic  promoters.  (A)  Natural  occurring  promoters  may 
have  hundreds  of  possible  TFBS.  In  the  ovarian  specific  promoter  1  (OSPI)  there  are  more  than  a  hundred 
possible  TFBS,  each  color  represents  a  different  TFBS  family.  (B)  Synthetic  promoter  for  Jun  (API)  contains 
several  TFBS  for  Jun  upstream  to  minimal  adenovirus  promoter  (lateADEp)  and  TSS.  (C)  the  synthetic  Jun(APl) 
we  constructed  is  activated  specifically  in  one  cell  line  (HCT)  and  not  the  others.  OSPI  appears  nonspecific, 
being  activated  in  most  cell  lines  shown.  Both  promoters  show  comparable  activation  levels  in  multiple  cancer 
(HCT,  MDA-453,  MDA-231,  OVCAR8)  and  normal  cell  lines  (MCF10A,  MCF12A,  adult  HDF). 
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Table  1:  Identification  of  promoters  that  are  active  in  prostate  cancers.  The  activity  of  14  breast 
cancer-specific  promoters  were  tested  on  2  prostate  cancer  cell  lines  (PC3,  DU145)  by  lentiviral 
transduction.  MFI  =  mean  fluorescence  intensity.  Fold  to  background  MFI  was  calculated  by  dividing 
the  MFI  of  individual  promoters  to  the  MFI  of  infection  control  cells.  2-3  promoters  exhibit  >  9  fold 
higher  activity  than  background  fluorescence  and  can  be  further  used  to  build  the  TIGRIS  system. 
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4.  IMPACT: 

What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

Our  data  utilizing  TIGRIS  system  for  treating  ovarian  cancer  had  demonstrated  the  feasibility  and  the 
potential  robust  efficacy  of  TIGRIS  system.  Furthermore,  the  sensor  library  and  screening  pipeline  we 
developed  enables  us  to  identify  sensors  for  various  types  of  tumor  cells.  If  the  adaptation  of  TIGRIS  system 
to  prostate  cancer  is  successful,  this  technology  would  represent  a  major  step  forward  in  our  ability  to 
eliminate  prostate  cancer,  especially  metastatic  and  recurrent  disease,  in  a  scalable  and  programmable 
fashion. 

What  was  the  impact  on  other  disciplines? 

Our  data  highlighted  the  potential  power  of  utilizing  gene  circuit  to  treat  unmet  clinical  needs,  such  as 
autoimmune  diseases. 

What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

Nothing  to  report. 

5.  CHANGES/PROBLEMS 

Changes  in  approach  and  reasons  for  change 

Nothing  to  report. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

Nothing  to  report. 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animals,  biohazards,  and/or  select 
agents 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  human  subjects 


Nothing  to  report. 
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Significant  changes  in  use  or  care  of  vertebrate  animals 

Nothing  to  report. 

Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  report. 

6.  PRODUCTS: 

Publications,  conference  papers,  and  presentations 
Journal  publications. 

Nissim  L,  Wu  MR  (Co-first  author),  Pery  E,  Binder-Nissim  A,  Suzuki  HI,  Stupp  D,  Wehrspaun  C, 
Tabach  Y,  Sharp  PA,  and  Lu  TK.  Synthetic  RNA-based  immunomodulatory  gene  circuits  for  cancer 
immunotherapy.  Cell  171  (5):  1 138- 1 150.el5,  2017;  acknowledgement  of  federal  support  (yes). 

Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  report. 

Other  publications,  conference  papers,  and  presentations. 

Nothing  to  report. 

Website(s)  or  other  Internet  site(s) 

Nothing  to  report. 

Technologies  or  techniques 

Nothing  to  report. 

Inventions,  patent  applications,  and/or  licenses 

Inventor  in  a  patent  based  on  my  postdoctoral  work,  submitted  by  Massachusetts  Institute  of  Technology: 

Lu,  Wu  MR  et  al.  (2016)  U.S.  TUMOR  IMMUNOTHERAPY.  62/181906,  M0656.70356 

Inventor  in  a  patent  based  on  my  postdoctoral  work,  submitted  by  Massachusetts  Institute  of  Technology: 

Lu,  Wu  MR  et  al.  (2017):  U.S.  SYNTHETIC  PROMOTERS.  62/470,754 


Other  Products 

Nothing  to  report. 
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7.  PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 


Name: 

Ming-Ru  Wu 

Project  Role: 

Postdoc  Researcher 

Researcher  Identifier 
(e.g.  ORCID  ID): 

orcid.org/0000-0002-2533-8484 

Nearest  person  month 
worked: 

12 

Contribution  to  Project: 

Dr.  Wu  has  been  performing  designing  and  optimizing 
gene  circuit  for  building  TIGRIS  system  for  prostate  cancer 

Funding  Support: 

Department  of  Defense  (W81XWH-16-1-0565). 
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SUMMARY 

Despite  its  success  in  several  clinical  trials,  cancer 
immunotherapy  remains  limited  by  the  rarity  of 
targetable  tumor-specific  antigens,  tumor-mediated 
immune  suppression,  and  toxicity  triggered  by  sys¬ 
temic  delivery  of  potent  immunomodulators.  Here, 
we  present  a  proof-of-concept  immunomodulatory 
gene  circuit  platform  that  enables  tumor-specific 
expression  of  immunostimulators,  which  could 
potentially  overcome  these  limitations.  Our  design 
comprised  de  novo  synthetic  cancer-specific  pro¬ 
moters  and,  to  enhance  specificity,  an  RNA-based 
AND  gate  that  generates  combinatorial  immunomod¬ 
ulatory  outputs  only  when  both  promoters  are  mutu¬ 
ally  active.  These  outputs  included  an  immunogenic 
cell-surface  protein,  a  cytokine,  a  chemokine,  and  a 
checkpoint  inhibitor  antibody.  The  circuits  triggered 
selective  T  cell-mediated  killing  of  cancer  cells,  but 
not  of  normal  cells,  in  vitro.  In  in  vivo  efficacy  assays, 
lentiviral  circuit  delivery  mediated  significant  tumor 
reduction  and  prolonged  mouse  survival.  Our  design 
could  be  adapted  to  drive  additional  immunomodu¬ 
lators,  sense  other  cancers,  and  potentially  treat 
other  diseases  that  require  precise  immunological 
programming. 

INTRODUCTION 

The  potential  of  cancer  immunotherapy  has  been  demonstrated 
in  several  clinical  trials,  yet  this  approach  is  still  limited  by  several 
major  challenges,  including  the  lack  of  targetable  tumor-spe¬ 


cific  antigens  and  tumor-mediated  immunosuppression.  For 
example,  chimeric  antigen  receptor  (CAR)-T  cells  and  bispecific 
antibodies  should  ideally  target  cell-surface  antigens  that  are 
exclusively  present  on  tumor  cells,  since  targeting  normal  tis¬ 
sues  can  result  in  severe  side  effects  (Morgan  et  al.,  201 0).  How¬ 
ever,  finding  highly  tumor-specific  surface  antigens  is  difficult, 
which  limits  the  range  of  targetable  tumors  (Klebanoff  et  al., 
2016).  In  addition,  even  when  ideal  targetable  tumor  antigens 
are  available,  tumor-mediated  immunosuppression  can  prevent 
successful  immunotherapy  by  disrupting  important  immunolog¬ 
ical  functions  that  are  necessary  for  effective  anti-tumor  immune 
responses  (Rabinovich  et  al.,  2007). 

Immunostimulatory  factors  can  be  utilized  to  overcome  these 
challenges,  but  off-target  activity  often  results  in  severe  toxicity. 
For  example,  surface  T  cell  engagers  (STEs)  are  artificial  immu¬ 
nogenic  cell-surface  proteins  that  bind  the  non-variable  regions 
of  the  T  cell  receptor  complex  (Liao  et  al.,  2000,  2003).  Conse¬ 
quently,  STE-expressing  cells  are  designated  forT  cell-mediated 
killing  regardless  of  T  cell  receptor  antigen  specificity  (see  Fig¬ 
ure  1 B  for  our  STE  design  based  on  membrane-anchored  anti- 
CD3  single-chain  variable  fragment  [scFv]).  Tumor-specific 
STE  expression  can  substitute  for  targetable  tumor  antigens 
but  must  be  constrained  to  cancer  cells  to  avoid  damage  to 
healthy  tissues  (Liao  et  al.,  2003).  Other  immunomodulators, 
such  as  chemokines,  cytokines,  and  immune  checkpoint  inhibi¬ 
tors,  can  assist  in  overcoming  tumor-mediated  immunosuppres¬ 
sion  but  have  caused  severe  side  effects  in  clinical  trials  (Lasek 
et  al.,  2014;  Leonard  et  al.,  1997).  Furthermore,  individual  immu¬ 
nomodulators  are  frequently  insufficient  on  their  own  to  yield 
strong  anti-tumor  efficacy  (Mahoney  et  al.,  2015).  Combinatorial 
immunotherapy  can  provide  significantly  stronger  efficacy  but 
also  increase  the  risk  and  severity  of  adverse  effects  (Boutros 
et  al.,  2016).  Tumor-localized  release  of  immunomodulators 
could  potentially  decrease  these  systemic  side  effects  and 
improve  the  therapeutic  efficacy,  but  is  challenging.  Thus, 
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Figure  1.  Immunomodulatory  Synthetic  Gene  Circuits  Are  Programmed  to  Be  Selectively  Activated  in  Cancer  Cells 

(A)  Panel  1 :  Immunomodulatory  synthetic  gene  circuits  are  designed  to  integrate  the  activity  of  two  tumor-specific  synthetic  promoters  (PI  and  P2)  with  an  RNA- 
based  AND  gate  mechanism  and  generate  combinatorial  immunomodulator  outputs  only  when  both  input  promoters  are  mutually  active.  When  activated,  the 
AND  gate  expresses  a  synthetic  transcription  factor  (GAD:  a  fusion  protein  consisting  of  the  GAL4  DNA  binding  domain  and  VP16  transcription  activating 
domain),  which  drives  the  co-expression  of  combinatorial  immunomodulators,  including  surface  T  cell  engagers  (STEs,  which  are  anti-CD3e  scFvs  displayed  on 
the  cell  surface)  and  secreted  CCL21 , 1  LI  2,  and  an  anti-PDI  antibody.  Panel  2:  The  circuits  are  triggered  to  express  immunomodulators  in  cancer  cells  (red),  but 
not  normal  cells  (blue).  Black-line  circuit  diagrams  indicate  circuits  are  active  and  gray-line  circuit  diagrams  indicate  circuits  are  inactive.  Panel  3:  Cancer-specific 
expression  of  combinatorial  immunomodulators  triggers  effective  T  cell-mediated  killing  of  the  cancer  cells. 

(B)  Schematic  drawing  of  STE  displayed  on  cell  surfaces. 


although  various  strategies  are  being  explored  to  overcome  this 
problem  (Neri  and  Sondel,  2016;  Shukla  and  Steinmetz,  2016; 
Tugues  et  al.,  2015),  further  development  is  required. 

Synthetic  biology  enables  the  creation  of  powerful  genetic 
biological  tools  for  studying,  diagnosing,  and  treating  disease 
(Bacchus  et  al.,  2013;  Kotula  et  al.,  2014;  Pardee  et  al.,  2016; 
Schukur  et  al.,  2015;  Xie  et  al.,  2016).  For  example,  synthetic 
tumor-targeting  gene  circuits  have  previously  been  designed 
to  target  tumor  cells  with  enhanced  specificity  but  have  not 
leveraged  the  anti-tumor  potential  of  the  immune  system,  which 
may  limit  their  efficacy  (Liu  et  al.,  201 4;  Morel  et  al.,  201 6;  Nissim 
and  Bar-Ziv,  2010;  Xie  et  al.,  2011).  We  hypothesized  that  gene 
circuits  could  be  designed  to  achieve  tumor-specific  production 
of  multiple  immunomodulators  and  thus  overcome  major  obsta¬ 
cles  in  current  immunotherapies  (Figure  1A). 

Here,  we  demonstrate  the  potential  of  gene  circuits  for  cancer 
immunotherapy.  Tumor-targeting  gene  circuits  often  require 
promoters  with  high  cell-state  specificity  to  function  as  cancer 
sensors,  which  can  be  hard  to  find  (Selvakumaran  et  al.,  2001; 


Yun  et  al.,  2008).  Moreover,  native  mammalian  promoters  can 
span  several  kb  and  are  therefore  difficult  to  fit  into  the  limited 
encoding  capacity  of  many  viral  vectors.  Thus,  we  constructed 
compact  synthetic  promoters  de  novo  that  exhibited  enhanced 
specificity  for  human  cancer  cells  compared  with  previously 
described  endogenous  cancer-specific  promoters. 

Discriminating  cancer  cells  based  on  the  integration  of  multi¬ 
ple  cancer  markers  significantly  enhances  targeting  specificity 
(Liu  et  al.,  2014;  Morel  et  al.,  2016;  Nissim  and  Bar-Ziv,  2010; 
Roybal  et  al.,  2016;  Xie  et  al.,  2011).  Thus,  to  enhance  the  circuit 
selectivity  over  a  single  synthetic  promoter,  we  designed  a  syn¬ 
thetic  RNA-based  circuit  that  integrates  the  activity  of  two  syn¬ 
thetic  ovarian-cancer-specific  promoters  and  generates  output 
proteins  only  when  both  promoters  are  decidedly  active,  analo¬ 
gous  to  a  Boolean  AND  gate.  Specifically,  input  promoter  1  (PI) 
regulates  an  auto-inhibitory  mRNAthat  encodes  the  output  pro¬ 
tein,  while  input  promoter  2  (P2)  regulates  an  RNA  molecule  that 
relieves  the  auto-inhibition  and  thus  enables  output  production. 
This  RNA-based  design  has  the  potential  to  reduce  the  off-target 
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immunogenic  signature  of  the  circuit  by  minimizing  the  expres¬ 
sion  of  foreign  proteins  in  normal  cells.  Immunomodulatory 
outputs  used  in  this  study  include  a  cell-surface  antigen  (STE), 
a  chemokine  that  promotes  T  cell  trafficking  to  tumor  sites 
(CCL21),  a  cytokine  that  enhances  T  cell  activation  and  function 
(IL12),  and  an  immune  checkpoint  inhibitor  (anti-PDI  antibody), 
hereafter  referred  to  as  the  SCIP  combination  (Lasek  et  al., 
2014;  Liao  et  al.,  2000, 2003;  Lin  et  al.,  201 4;  Postow  et  al.,  201 5). 

We  characterized  our  circuits  with  in  vitro  specificity  assays 
using  normal  and  cancer  cell  lines  and  in  vivo  efficacy  assays  us¬ 
ing  a  disseminated  ovarian  cancer  model,  which  recapitulates 
advanced  human  ovarian  cancer  (Lengyel  et  al.,  2014).  In  addi¬ 
tion,  we  demonstrated  that  our  circuit  can  be  readily  modified 
to  target  breast  cancer  cells  by  replacing  the  synthetic  input  pro¬ 
moters.  Our  circuit  design  is  modular  and  enables  facile  con¬ 
struction,  testing,  and  optimization  of  arbitrary  promoter  inputs, 
circuit  components,  and  outputs.  The  outputs  can  be  any  genet¬ 
ically  encodable  element,  including  intracellular,  secreted,  and 
cell-surface  proteins.  This  flexibility  will  enable  circuits  to  be 
optimized  using  advanced  tumor  models  and  patient-derived 
tumors,  which  will  be  important  for  translation  into  clinical  use. 
Thus,  our  programmable  immunomodulatory  gene  circuits 
have  broad  potential  utility  for  cancer  immunotherapy  and  study¬ 
ing  tumor  immunology. 

RESULTS 

RNA-Based  AND  Gate  Circuit 

We  first  designed  and  optimized  the  RNA-based  AND  gate.  We 
engineered  two  genetic  modules  such  that  each  was  regulated 
by  a  separate  promoter  (Figures  2A  and  2B).  Module  1  consists 
of  PI ,  which  regulates  the  expression  of  an  RNA  that  encodes 
two  exons  of  an  output  protein  mKate2  (mK-ExI,  mK-Ex2), 
separated  by  a  synthetic  intronic  microRNA  (miRNA)  (miRI),  as 
well  as  perfect  match  miRI  binding  sites  (BS[Pe])  downstream 
to  the  second  exon.  Following  splicing,  the  exons  are  assembled 
into  a  mature  output  mRNA  that  contains  BS(Pe)s  in  the  3'  un¬ 
translated  region  (UTR),  while  the  intron  is  processed  by  the 
cell  to  produce  miRI.  Consequently,  when  only  PI  is  active, 
miRI  binds  the  BS(Pe)  sites  in  the  output  mRNA  and  targets  it 
for  degradation.  Therefore,  module  1  constitutes  an  auto-inhibi- 
tory  loop  for  output  protein  production  (Figure  2A).  Module  2 
consists  of  P2,  which  regulates  a  miRNA  sponge  that  consists 
of  multiple  bulged  complementary  miRI  binding  sites  (BS[B]) 
and  thus  acts  as  a  miRI  shunt  that  relieves  the  auto-inhibitory 
nature  of  module  1  when  it  is  expressed  (Figures  2A  and  2B). 
In  this  system,  the  output  protein  is  expressed  at  a  high  level 
only  when  the  promoters  regulating  both  modules  are  mutually 
active.  However,  when  only  one  of  the  promoters  is  active  or 
when  none  of  the  promoters  are  active,  output  protein  expres¬ 
sion  is  low. 

We  systematically  optimized  various  aspects  of  the  miRNA- 
based  mechanism  (Figures  3,  SI ,  S2,  S3,  and  STAR  Methods) 
(Auyeung  et  al.,  2013;  Ebert  and  Sharp,  2010;  Fellmann  et  al., 
2013;  Suzuki  et  al.,  2015).  In  module  2,  we  characterized 
how  the  miRNA  binding  site  (BS)  sequence,  expression  level, 
and  sponge  architecture  determines  the  efficiency  of  the  miRI 
sponge  (Figures  3,  SI ,  and  S2).  We  defined  four  different  states 


for  this  circuit:  in  state  [0,0],  both  module  1  and  module  2  are 
inactive;  in  state  [0,1],  module  1  is  inactive,  while  module  2  is 
active;  in  state  [1,0],  module  1  is  active,  while  module  2  is  inac¬ 
tive;  and  in  state  [1,1],  both  module  1  and  module  2  are  active 
(Figure  2B).  The  ON-OFF  ratio  for  the  AND  gate  circuit  was 
defined  as  the  mKate2  output  level  in  the  presence  versus  the 
absence  of  the  sponge  (state  [1 ,1]  versus  state  [1 ,0]  in  Figure  2B). 
Ultimately,  we  chose  Sponge  7:  H2A1p-EC1-10BZRANB1-EC2- 
5B  (Figure  S2A)  as  the  final  “optimized  sponge”  design  because 
it  yielded  the  highest  ON-OFF  ratio  for  the  AND  gate  circuit 
(-6-fold,  Figure  S2A). 

In  module  1,  we  optimized  the  synthetic  intronic  pri-miRNA 
backbone  in  the  self-inhibitory  mKate2  transcript  for  enhanced 
miRI  production  efficiency  (Greber  et  al.,  2008;  Nissim  et  al., 

2014;  Xie  et  al.,  2011).  The  absolute  mKate2  signal  in  the  OFF 
(no  sponge,  state  [1 ,0])  and  ON  states  (optimized  sponge,  state 
[1,1])  were  both  affected  by  the  pri-miRNA  designs  we  tested 
(Mv1-Mv3,  Figure  S3).  Mvl  exhibited  an  OFF-ON  shift  from 
—1 400  to  —8400  mean  fluorescence  intensity  (MFI)  (5.8-fold  acti¬ 
vation  ratio),  Mv2  from  —170  to  —500  MFI  (3-fold  activation),  and 
Mv3  from  —380  to  —2600  MFI  (6.8-fold  activation)  (Figures  S3D 
and  S3E).  We  selected  the  Mv3  architecture  for  further  experi¬ 
mentation  because  it  had  low  mKate2  level  in  the  OFF  state,  sub¬ 
stantial  absolute  mKate2  levels  in  the  ON  state,  and  the  highest 
fold  activation.  However,  we  note  that  different  designs  could  be 
more  suitable  for  other  targeting  scenarios.  For  example,  when 
using  effector  outputs  with  low  potency,  Mvl  could  be  chosen 
to  achieve  high  absolute  levels  of  expression  in  the  ON  state. 
Conversely,  Mv2  could  be  suitable  for  expressing  extremely 
potent  outputs  for  which  even  low  background  output  activity 
could  be  detrimental. 

Synthetic  Cancer-Specific  Promoters 

To  generate  synthetic  ovarian-cancer-specific  promoters,  we 
identified  transcription  factors  (TFs)  that  are  overexpressed  in 
ovarian  cancer  cells  by  comparing  gene  expression  in  primary 
human  ovarian  tumors  versus  healthy  controls,  and  the  respec¬ 
tive  binding  motifs  of  these  TFs,  using  publically  available  data¬ 
bases  (STAR  Methods).  In  each  synthetic  promoter  S(TF)p,  mul¬ 
tiple  binding  motifs  for  a  single  cancer-specific  TF  were  encoded 
in  tandem  upstream  to  a  synthetic  minimal  promoter  derived 
from  the  major  adenoviral  late  promoter  (Figure  4A).  This  design 
approach  resulted  in  compact  promoters  with  evenly  distributed 
and  dense  TF-BSs  corresponding  to  a  small  number  of  TF  fam¬ 
ilies.  For  example,  analysis  of  the  1  -kb  upstream  of  the  transcrip¬ 
tional  start  site  (TSS)  of  the  native  tumor-specific  promoter 
H2A1  p  (Rogakou  et  al.,  1998)  revealed  -82  possible  BSs  for 
40  different  TF  families,  with  a  mean  of  —8  TF-BSs  per  100  bp. 
In  comparison,  the  synthetic  S(E2F1)p  promoter,  comprised  of 
E2F1  BSs,  had  only  44  possible  BSs  for  7  different  TF  families, 
with  a  mean  of  —16  TF-BSs  per  100  bp  (Table  SI  and  STAR 
Methods). 

We  assayed  the  specificity  and  activity  of  natural  and  artificial 
promoters  in  human  OVCAR8  (OV8)  ovarian  cancer  cells  and 
normal  primary  cells,  including  primary  adult  human  dermal  fi¬ 
broblasts  (aHDF)  and  primary  human  ovarian  epithelium  (HO V). 
Several  synthetic  promoters,  including  S(E2F1)p  and  S(cMYC) 
p  (comprised  of  E2F1  and  cMYC  BSs,  respectively),  were  highly 
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Figure  2.  The  Basic  RNA-Only  Single-Output  AND  Gate  Design 

(A)  Module  1  of  the  AND  gate  is  designed  as  an  auto-inhibitory  loop  such  that  it  represses  its  own  output.  Module  2  is  designed  to  inhibit  the  auto-inhibition  of 
module  1 .  Module  1  and  module  2  are  regulated  by  cancer-specific  promoters  PI  and  P2,  respectively.  The  output  from  module  1  is  expressed  at  a  high  level  only 
when  both  PI  and  P2  are  active,  which  enhances  the  tumor  specificity  of  the  circuit. 

(B)  All  four  possible  input  states  and  their  respective  output  states  are  shown  for  the  RNA-only  single-output  AND  gate.  Input  states  are  defined  within  the  square 
brackets  by  whether  module  1  and  module  2  are  active,  where  0  means  inactive  and  1  means  active.  In  state  [1,0],  PI  is  active  and  the  mKate2  transcript  is 
expressed  (mK-ExI  and  mK-Ex2  denote  mKate2  exonl  and  mKate2  exon2,  respectively).  However,  this  transcript  encodes  miRI  within  the  mKate2  gene  that 
inhibits  the  mKate2  transcript  by  targeting  BS(Pe)s  in  the  3'  UTR.  Thus,  the  mKate2  levels  are  minimal.  In  state  [0,1],  PI  is  inactive,  so  the  output  protein  mKate2  is 
not  expressed.  In  state  [0,0],  neither  PI  nor  P2  are  active,  and  thus,  the  output  protein  mKate2  is  not  expressed.  In  state  [1,1],  the  P2  promoter  expresses  a  sponge 
for  miRI  that  is  based  on  BS(B)s.  This  enables  sequestration  of  miRI  away  from  inhibiting  the  mKate2  transcript  expressed  by  PI,  thus  allowing  for  mKate2 
expression. 


active  in  OV8  ovarian  cancer  cells,  but  not  in  normal  cells. 
These  synthetic  promoters  exhibited  greater  activation  ratios 
compared  with  known  native  tumor-specific  promoters  such 
as  SSXIp  (Gure  et  al.,  1997)  and  H2A1p  (Rogakou  et  al.,  1998) 

(Figure  4B  for  selected  promoters). 


For  example,  SSXIp  exhibited  9-fold  activation  in  OV8 
compared  to  aHDF  and  17-fold  activation  compared  to  HOV 
(Figure  4B).  H2A1p  had  2-fold  activation  in  OV8  compared  to 
aHDF  and  5.6-fold  compared  to  HOV.  The  synthetic  promoter 
S(E2F1)p  had  224-fold  induction  in  OV8  versus  aHDF  and 
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Figure  3.  Flowchart  of  Circuit  Optimization  Strategies  for  the  RNA- 
Based  AND  Gate 

Multiple  versions  of  the  miRI  production  backbone  in  module  1 ,  as  well  as  the 
miRI  BS  sequences  and  sponge  architecture  in  module  2,  were  systematically 
tested.  The  optimized  circuit  utilized  the  miR-E  backbone  (Mv3,  Table  S2)  in 
module  1  to  generate  miRI  and  a  sponge  architecture  “Sponge  7”  in  module  2 
that  encoded  ten  optimized  BS(B)s  within  an  enhanced  cyan  fluorescent 
protein  (ECFP)  gene  (EC1-10BZRAnbi-EC2)  and  five  additional  BS(B)s  in  the 
downstream  3'  UTR  (Figure  S2A).  The  mKate2  output  levels  for  the  optimized 
circuit  with  a  control  sponge  with  no  miRI  BSs  (“Mv3  +  Cont.”,  representing 
the  [1,0]  state  in  figure  2B),  and  the  optimized  circuit  with  the  optimized 
sponge  (“Mv3  +  Sponge”,  representing  the  [1 ,1]  state  in  Figure  2B)  are  shown. 
Error  bars  represent  SEM,  n  =  3  biological  replicates.  TuD  =  tough-decoy 
architecture. 

See  also  Figures  SI ,  S2,  and  S3;  Table  S2;  and  STAR  Methods. 


570-fold  activation  in  OV8  compared  to  HOV,  with  absolute 
maximal  expression  levels  comparable  to  the  native  promoters. 
S(cMYC)p  exhibited  15-fold  induction  in  OV8  versus  aHDF  and 
36-fold  activation  in  OV8  compared  to  HOV.  Since  the  maximal 
activity  of  S(cMYC)p  was  lower  than  that  of  S(E2F1)p  in  OV8 
cells,  we  selected  S(cMYC)p  to  drive  expression  of  module  1 
and  S(E2F1)p  to  drive  the  expression  of  module  2.  This  was 
done  so  that  the  sponge  could  be  expressed  at  a  greater  level 
than  the  miRI -encoding  construct,  which  is  important  for 
achieving  a  high  activation  ratio  with  our  circuit  (Figure  S2A). 

Multi-output  AND  Gate  for  Combinatorial 
Immunomodulation 

To  facilitate  the  expression  of  combinatorial  immunomodulators 
from  our  circuit,  we  replaced  mKate2  in  module  1  with  a  synthetic 
transcription  factor  (GAD)  that  binds  to  its  cognate  synthetic 
promoter  and  activates  the  expression  of  downstream  genes 
(Figure  5A).  This  GAD  was  previously  made  by  fusing  the  yeast 


GAL4  DNA-binding  domain  to  the  viral  VP16  transcriptional 
activation  domain.  Specifically,  we  used  the  synthetic  promoter 
S(cMYC)p  to  express  module  1  (S(cMYC)p-GADEXr[miR1]cons- 
GADEX2-11Pe)-  Module  1  encoded  a  self-inhibitory  RNA  tran¬ 
script  with  miRI  as  an  intron  within  the  GAD  coding  sequence 
followed  by  11  tandem  BS(Pe)s  in  the  3'  UTR.  In  the  module  2 
vector,  the  synthetic  promoter  S(E2F1)p  regulated  the  expres¬ 
sion  of  the  optimized  sponge  architecture  ( Sponge  7:  EC1- 
1  0BZranbi-EC2-5B).  Finally,  we  encoded  an  mKate2  output  un¬ 
der  the  regulation  of  an  artificial  promoter  targeted  by  GAD  that 
contains  GAL4-BSs  upstream  to  a  synthetic  minimal  late  adeno¬ 
viral  promoter  (GALp),  namely  module  3  (Figure  5A).  We  hypoth¬ 
esized  that  increasing  the  number  of  GAL4-BSs  in  the  GALp 
would  increase  output  expression  levels.  We  therefore  built 
3  different  versions  of  GALp,  either  with  5,  8,  or  14  tandem 
GAL4-BSs  (G5p,  G8p,  and  G14p,  respectively).  We  further  sur¬ 
mised  that  output  levels  in  the  OFF  state  could  be  lowered  by 
adding  miRI-BS(Pe)  to  the  mKate2  output,  thus  allowing  miRI 
to  reduce  both  GAD  and  mKate2  levels  directly  when  the  sponge 
is  not  active  (Figure  5A).  Thus,  we  built  three  additional  versions 
of  the  output  vector  in  which  a  single  BS(Pe)  was  encoded  in  the 
3'  UTR  of  the  output  mKate2  transcript  (G5-Pe,  G8-Pe  and  G14- 
Pe,  respectively).  This  design  architecture  enables  convenient 
tuning  of  each  output’s  expression  levels  based  on  the  promoter 
and  whether  one  Pe  was  incorporated  into  the  output  transcript. 

We  first  characterized  the  multi-output  AND-gate  architec¬ 
tures  by  modeling  scenarios  in  which  neither  (state  [0,0]),  only 
one  (states  [0,1]  or  [1 ,0]),  or  both  (state  [1,1])  cancer-specific  pro¬ 
moters  are  active  in  OV8  cells  (Figure  5B  and  STAR  Methods). 
The  G8-Pe  circuit  exhibited  negligible  mKate2  levels  in  states 
[0,0]  and  [0,1]  and  high  8-fold  mKate2  activation  in  state  [1,1] 
over  the  minimal  mKate2  expression  in  state  [1,0]  and  was 
thus  chosen  for  further  use. 

We  next  compared  mKate2  expression  from  the  G8-Pe  circuit 
delivered  via  lentivirus  in  OV8  cells  versus  normal  cells,  including 
primary  aHDF,  HOV,  and  human  T  cells.  The  circuit  exhibited 
higher  mKate2  expression  in  OV8  versus  aHDF  and  HOV  in  the 
state  [1,1]  circuit  configuration  (4-fold  and  71 -fold,  respectively; 
Figure  5C).  We  next  determined  the  activity  of  the  S(cMYC)p  and 
S(E2F1)p  promoters  individually,  as  well  as  the  activity  of  the 
G8-Pe  circuit  with  S(cMYC)p  and  S(E2F1)p  promoters  as  inputs, 
in  primary  human  T  cells.  We  observed  that  S(E2F1)p  was 
significantly  activated  in  primary  human  T  cells.  However,  since 
S(cMYC)p  was  not  active  in  these  cells,  the  integrated  output 
generated  by  the  G8-Pe  circuit  was  30-fold  lower  compared  to 
the  activity  of  S(E2F1)p  on  its  own  (Figures  5D  and  S4A).  Further¬ 
more,  we  tested  the  specificity  of  this  design  using  additional 
cell  lines,  including  non-tumorigenic  human  ovarian  fibro¬ 
blasts  (HOFs),  human  ovarian  microvasculature  epithelial  cells 
(HOMECs),  colonic  epithelial  cells  (CCD-841 -CoN),  and  mam¬ 
mary  epithelial  cells  (MCF-1 0A  and  MCF-1 2A),  as  well  as  induc¬ 
ible  pluripotent  stem  cells  (iPSCs).  The  circuit  achieved  high 
specificity  for  OV8  cells,  with  10-fold  to  800-fold  higher  circuit 
output  in  OV8  cells  versus  these  additional  cell  lines  (Figure  S4B). 
These  data  demonstrate  the  stringency  of  our  circuit  and  the  po¬ 
tential  of  our  circuit  design  to  enhance  therapeutic  specificity  by 
reducing  the  probability  of  false-positive  targeting  compared  to 
targeting  with  a  single  promoter. 
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Figure  4.  Synthetic  Tumor-Specific  Pro¬ 
moters  Exhibited  High  Tumor  Specificity 

(A)  Design  of  synthetic  tumor-specific  promoters, 
S(TF)p.  Multiple  BSs  for  a  TF  that  is  overexpressed 
in  cancer  (TF-BS)  were  cloned  in  tandem  upstream 
of  a  minimal  late  adenovirus  promoter  (lateADEp). 

(B)  Synthetic  tumor-specific  promoters  regulating 
mKate2  exhibited  enhanced  tumor  specificity. 
Representative  native  tumor-specific  promoters, 
SSXIp  and  FI2A1p,  are  shown  for  comparison. 
Parentheses  indicate  the  TF-BSs  used  to  construct 
each  synthetic  tumor-specific  promoter.  aFIDF  and 
FIOV  are  normal  primary  cells.  OV8  is  a  human 
ovarian  cancer  cell  line.  Error  bars  represent  SEM, 
n  =  3  biological  replicates. 

See  also  Table  SI. 


Finally,  to  demonstrate  the  modularity  of  our  design,  we  exam¬ 
ined  whether  the  G8-Pe  circuit  could  be  adapted  to  target  breast 
cancer  cells  by  simple  modification  of  the  input  promoters 
(Figure  5E).  For  this  purpose,  we  created  two  additional  synthetic 
promoters,  S(USF1)p  and  S(MAFK)p,  that  were  designed  to  spe¬ 
cifically  target  the  MDA-MB-453  breast  cancer  cell  line,  but  not 
the  MCF-10A  non-tumorigenic  mammary  epithelial  cell  line. 
We  compared  the  mKate2  outputs  generated  by  each  promoter 
individually  versus  the  G8-Pe  circuit  in  both  cell  lines.  This  new 
design  generated  a  high  output  only  in  MDA-MB-453  cells, 
demonstrating  the  versatility  of  our  circuit  architecture  and  its 
ability  to  be  adapted  to  a  variety  of  different  therapeutic  sce¬ 
narios  and  applications. 

STEs  Trigger  T  Cell-Mediated  Cancer-Cell  Killing  and 
IFN-y  Secretion 

Next,  we  sought  to  encode  immunotherapies  as  the  outputs  of 
our  synthetic  cancer-specific  circuits.  An  STE  was  constructed 
by  fusing  an  scFv  targeting  human  CD3e  with  an  inert  membrane 
anchor  derived  from  the  Duffy  antigen/receptor  for  chemokines 
(DARC)  of  erythrocytes  (von  Nickisch-Rosenegk  et  al.,  2012) 
(Figure  IB).  We  replaced  mKate2  with  STE  as  the  output  of  the 
G8-Pe  circuit  and  validated  the  therapeutic  specificity  of  this  cir¬ 
cuit  in  vitro.  We  observed  20-  to  1 1 8-fold  higher  STE  expression 
on  circuit-transduced  OV8  cells  than  on  aHDF  or  HOV  cells  (Fig¬ 
ure  6A).  We  also  observed  robust  T  cell-mediated  killing  of  cir¬ 
cuit-transduced  OV8  cells,  which  was  significantly  higher  than 
in  aFIDF  or  FIOV  cells  (Figures  6B  and  S4C-S4E).  There  was  min¬ 
imal  T  cell-mediated  killing  of  non-STE-expressing  cells  (Figures 
S4C-S4E),  which  could  have  been  caused  by  our  experimental 
conditions,  including  major  histocompatibility  complex  (MHC) 
mismatch  between  T  cells  and  target  cells,  as  well  as  T  cell- 
mediated  killing  with  non-STE-based  mechanisms.  In  addition 
to  direct  tumor  killing,  robust  interferon  (IFN)-y  secretion  by 
T  cells  is  critical  for  strong  therapeutic  efficacy.  We  observed 


strong  and  specific  IFN-y  secretion  by 
T  cells  only  when  they  were  co-cultured 
with  circuit-transduced  OV8  cells  (Figures 
6C  and  S4F),  which  was  20-  to  28-fold 
higher  than  when  T  cells  were  co-cultured 
with  circuit-transduced  aHDF  or  HOV 
cells.  We  also  observed  minimal  OV8  killing  and  IFN-y  secretion 
by  T  cells  when  module  2  was  removed  from  the  G8-Pe  circuit 
(state  [1,0])  to  model  the  potential  scenario  in  which  one  of  the 
promoters  exhibits  some  off-target  activity  (Figures  6A-6C  and 
S4C-S4F). 

We  also  performed  a  control  experiment  in  which  STE  expres¬ 
sion  was  regulated  by  the  human  ubiquitin  C  (hubC)  promoter, 
which  is  constitutively  active  in  OV8,  aHDF,  and  HOV  cells.  We 
examined  the  sensitivity  of  these  cells  to  STE-mediated  killing 
and  the  cells’  ability  to  trigger  IFN-y  production  by  T  cells  (Fig¬ 
ures  S4G-S4I).  We  found  that  while  STE  expression  levels  in 
OV8  cells  were  ~3-fold  higher  than  in  aHDF  and  HOV  cells, 
STE-mediated  killing  was  comparable  in  aHDF  and  OV8  and 
~2-fold  higher  in  HOV.  In  addition,  all  three  STE-expressing 
cell  lines  triggered  comparably  strong  IFN-y  production  by 
T  cells.  These  results  highlight  the  stringency  of  our  AND  gate, 
which  mediates  significantly  higher  OV8  specificity  over  these 
other  cell  lines. 

We  next  examined  whether  STE  expression  could  trigger 
T  cell-mediated  anti-tumor  activity  in  vivo.  OV8  cells  were  trans¬ 
duced  with  lentiviruses  encoding  a  doxcycline-inducible  STE 
expression  system  and  then  injected  intraperitoneally  (i.p.)  into 
immunodeficient  NOD.Cg-Pr/cc/cSGid  ll2rgtm1Wi'/ SzJ  (NSG)  mice 
on  day  0.  Since  NSG  mice  lack  T  cells,  human  T  cells  were 
also  periodically  injected  i.p.  Specifically,  STE  expression  was 
induced  at  day  14,  followed  by  human  T  cell  injections  at  days 
1 6  and  31 .  In  the  absence  of  T  cells,  tumor  growth  was  similar  be¬ 
tween  the  STE-induced  and  non-induced  groups.  STE  expres¬ 
sion  led  to  a  6-fold  decrease  in  tumor  burden  in  mice  treated 
with  T  cells  compared  to  non-treated  mice  (Figures  S4J-S4L), 
demonstrating  that  the  anti-tumor  activity  of  STEs  was  depen¬ 
dent  on  T  cells.  However,  consistent  with  previous  studies 
(Liao  et  al.,  2003),  STE  expression  on  its  own  was  unable  to 
completely  abrogate  tumor  growth,  even  when  expressed  in 
high  levels  in  nearly  all  the  tumor  cells  (Figure  S4L). 
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Figure  5.  Multi-output  AND  Gate  Architecture  and  Performance 

(A)  The  incorporation  of  a  synthetic  transcription  factor  (GAD)  as  the  output  of  module  1  means  that  it  is  only  expressed  at  high  levels  when  both  PI  and  P2  are  high 
(state  [1 ,1]).  GAD  enables  facile  tuning  of  downstream  gene  expression  (mKate2  •  •  •  Output  N). 

(B)  mKate2  expression  can  be  modulated  by  modifying  the  number  of  GAD  BSs  in  the  synthetic  GALp  expressing  mKate2  (G5,  G8,  G14  for  5,  8,  14  GAD  BSs  in 
GALp)  and  the  presence  of  miRNA  BSs  at  the  mKate2  3'  UTR  (G5-Pe,  G8-Pe,  G14-Pe  indicate  transcripts  with  a  BS[Pe]).  Note  the  circuit  diagram  depicts  the 
version  with  a  perfect  miRI  BS  in  the  mKate2  transcript.  State  [0,0]  indicates  cells  containing  the  negative  control  modules,  module  Icon  and  module  2con  (see 
STAR  Methods  for  details).  State  [0,1]  indicates  cells  containing  module  Icon  and  module  2.  State  [1,0]  indicates  cells  containing  module  1  and  module  2con. 
State  [1,1]  indicates  cells  containing  module  1  and  module  2.  All  cell  states  tested  also  contained  their  respective  module  3. 

(C)  The  G8-Pe  AND  gate  architecture  triggered  high  mKate2  expression  in  human  ovarian  cancer  cells  (OV8)  but  not  in  normal  cells  in  state  [1 ,1]. 

(D)  Our  circuit  design  prevents  potential  off-target  effects  in  primary  human  T  cells.  The  S(E2F1)p  promoter  was  active  in  T  cells  but  S(cMYC)p  was  inactive  in 
T  cells.  The  G8-Pe  circuit  triggered  minimal  mKate2  output  in  T  cells,  thus  minimizing  the  potential  off-target  effects  that  would  have  been  observed  if  the  S(E2F1)p 
was  used  to  directly  drive  mKate2  expression. 


(legend  continued  on  next  page) 
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Tumor-Specific  Expression  of  Multiple 
Immunotherapeutic  Outputs  Triggers  Effective  T  Cell- 
Dependent  Tumor  Killing 

Since  combination  therapy  has  shown  enhanced  therapeutic  ef¬ 
ficacy  over  monotherapies  (Lasek  et  al.,  2014;  Mahoney  et  al., 

2015;  Tugues  et  al.,  2015)  and  STE  expression  alone  was  insuf¬ 
ficient  to  achieve  strong  anti-tumor  effects  in  vivo,  we  engineered 
our  optimized  cancer-specific  gene  circuit  to  generate  multiple 
immunomodulators  as  outputs.  To  demonstrate  the  capability 
of  the  circuit  to  generate  various  types  of  commonly  used  immu¬ 
nomodulators,  we  implemented  the  SCIP  combination. 

First,  we  transduced  OV8  cells  with  lentiviruses  encoding  the 
G8-Pe  AND  gate  and  the  SCIP  combination.  We  validated  the 
selective  co-production  of  STE,  CCL21,  IL12,  and  an  anti-PDI 
antibody  (Ab)  from  OV8  cells,  but  not  from  aHDF  or  HOV  cells, 
in  vitro  (Figure  6D-6G).  Next,  we  validated  the  functionality  of 
our  SCIP  circuit  with  an  in  vivo  model  of  human  ovarian  cancer. 
We  engineered  the  G8-Pe  circuit  to  produce  a  control  output, 
STE  only,  or  the  SCIP  combination  in  OV8  cells.  We  then  estab¬ 
lished  i.p.  disseminated  ovarian  cancers  by  injecting  these  OV8 
cells  into  the  peritoneal  cavity  of  NSG  mice.  This  experimental 
setup  was  used  to  model  the  functionality  of  the  circuit 
throughout  the  course  of  the  treatment  and  its  effects  on  circuit 
efficacy.  For  example,  the  activity  of  our  input  promoters  could 
be  altered  in  response  to  environmental  signals  or  due  to  tumor 
heterogeneity,  resulting  in  circuit  inactivation  and  reduced 
efficacy. 

In  the  absence  of  T  cells,  all  tumors  grew  robustly  (Figure  S5A). 
When  administered  with  T  cells  i.p.,  tumor  sizes  in  the  control 
output  and  STE-only  groups  were  significantly  larger  than  the 
STE+CCL21 +IL1 2+anti-PD1  Ab  (SCIP)  combination,  in  which  tu¬ 
mors  were  undetectable  in  most  mice  (Figures  7A,  S5A,  and 
S5B).  Significantly  prolonged  mice  survival  was  also  observed 
in  the  SCIP  group  (Figure  7B).  Thus,  combinatorial  immunother¬ 
apies  expressed  from  within  tumor  cells  can  mediate  strong  anti¬ 
tumor  effects. 

Next,  we  sought  to  confirm  that  the  SCIP-expressing  circuit 
could  mediate  robust  efficacy  even  when  it  was  only  introduced 
into  a  small  fraction  of  tumor  cells.  This  experiment  was  used  to 
model  a  common  problem  with  tumor-targeting  therapies  in 
which  the  therapeutic  circuit  is  activated  only  in  a  fraction  of 
the  tumor.  For  example,  this  phenomenon  can  be  mediated  by 
modulation  of  gene  expression,  cellular  heterogeneity  within  tu¬ 
mors,  or  partial  circuit  delivery.  We  mixed  wild-type  OV8  cells 
with  OV8  cells  containing  the  SCIP-expressing  circuit  in  defined 
ratios.  Even  when  only  15%  or  30%  of  the  overall  tumor  was 
composed  of  SCIP-expressing  OV8  cells,  significant  tumor 
reduction  and  survival  was  still  achieved,  but  only  in  the  pres¬ 
ence  of  T  cells  (Figures  7C,  7D,  S6A,  and  S6B).  These  data  sug¬ 
gest  that  immunomodulatory  gene  circuits  may  be  able  to 
mediate  strong  anti-tumor  effects  despite  incomplete  delivery 
in  future  applications. 


Finally,  we  directly  tested  the  efficacy  of  our  SCIP  circuit  within 
in  vivo  delivery  settings.  To  minimize  the  total  number  of  viruses 
required  to  encode  the  entire  SCIP  circuit  and  thus  simplify  the 
delivery  process,  we  combined  the  CCL21  and  IL1 2  output  mod¬ 
ules  into  a  single  lentiviral  construct  by  encoding  both  outputs, 
separated  by  a  P2A  peptide  (Szymczak  et  al.,  2004),  under  the 
regulation  of  a  single  G8  output  promoter.  We  similarly  combined 
the  anti-PDI  Ab  heavy-chain  (HC)  and  light-chain  (LC)  modules 
into  a  single  lentivirus.  Mice  with  i.p.  disseminated  OV8  tumors 
were  injected  i.p.  with  either  a  control  lentivirus  mix  that  con¬ 
tained  module  1 ,  module  2,  and  a  negative  control  output  or  a 
therapeutic  virus  mix  that  contained  module  1,  module  2,  and 
the  SCIP  outputs.  Human  T  cells  were  also  periodically  injected 
i.p.  Under  these  settings,  we  observed  significantly  reduced  tu¬ 
mor  burden  and  enhanced  survival  in  mice  treated  with  the  SCIP 
circuit  versus  the  control  unit  (Figures  7E,  7F,  S7A,  and  S7B). 
These  data  indicate  that  in  vivo  delivery  of  our  immunomodula¬ 
tory  gene  circuits  has  the  potential  to  trigger  anti-tumor  effects, 
even  when  the  components  of  our  gene  circuit  are  encoded  on 
multiple  vectors. 

DISCUSSION 

In  this  study,  we  present  a  proof-of-concept  demonstration  that 
versatile  genetic  circuits  can  be  engineered  for  highly  specific 
production  of  combinatorial  immunomodulators  from  within  can¬ 
cer  cells.  These  circuits  were  delivered  in  vivo  to  trigger  robust 
anti-tumor  responses  and  increased  survival  of  treated  mice. 
Our  modular  design  strategy  enables  the  facile  optimization 
required  for  future  clinical  translation  and  targeting  of  additional 
tumor  types. 

However,  several  issues  must  be  addressed  to  further 
advance  this  approach  toward  therapeutic  applications.  Design 
and  testing  of  circuits  will  likely  need  to  be  performed  in  addi¬ 
tional  models,  such  as  patient-derived  tumor  cells  and  a  panel 
of  normal  tissues.  The  potential  off-target  activity  of  circuit  com¬ 
ponents,  such  as  miRI  and  the  miRNA  sponge,  should  also 
be  examined,  for  example,  via  transcriptomics.  In  vitro  assays 
only  allow  for  testing  of  short-term  specificity,  and  thus,  future  ef¬ 
forts  should  also  examine  long-term  specificity  and  side  effects 
within  in  vivo  settings. 

Furthermore,  the  NSG  mice  supplemented  with  human  T  cells 
that  we  utilized  in  this  study  have  several  drawbacks  that  limit  our 
ability  to  fully  evaluate  safety  and  efficacy.  In  this  model,  tumor 
cells  are  of  human  origin,  but  normal  tissues  are  of  mouse  origin, 
thus  making  it  challenging  to  characterize  targeting  specificity 
for  human  applications.  Second,  NSG  mice  injected  with  human 
T  cells  eventually  develop  life-threatening  graft-versus-host  dis¬ 
ease  (Covassin  et  al.,  2011),  making  it  difficult  to  evaluate  long¬ 
term  safety  concerns.  Third,  these  mice  only  contain  a  partial  hu¬ 
man  immune  system,  established  by  human  T  cell  injections,  and 
may  not  be  sufficient  to  determine  the  exact  optimal  combination 


(E)  Our  circuit  design  can  be  readily  adapted  to  distinguish  a  different  tumor  type  from  its  normal  counterparts.  By  using  two  different  promoters,  the  G8-Pe  circuit 
triggered  a  strong  output  in  the  breast  cancer  cell  line  M DA- MB -453  but  minimal  output  in  the  non-tumorigenic  breast  epithelial  cell  line  MCF-10A. 

Error  bars  represent  SEM,  n  =  3  biological  replicates  for  all  experiments,  except  n  =  6  biological  replicates  for  OV8  group  in  (C). 

See  also  Figure  S4A. 
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of  immunomodulatory  outputs  required  for  efficacy  in  humans. 
Therefore,  future  work  will  be  needed  to  optimize  and  charac¬ 
terize  specific  immunomodulatory  combinations  in  additional 
animal  models,  including  humanized  mice,  syngeneic  mouse 
models  with  intact  immune  systems,  and  patient-derived  xeno¬ 
graft  models,  which  can  account  for  tumor  heterogeneity. 

Improved  delivery  strategies  will  be  needed  to  convey  our 
immunomodulatory  gene  circuits  into  tumors,  such  as  non- 
viral  carriers  or  oncolytic  viruses.  Yet,  given  that  only  15%  of 
tumor  cells  had  to  express  the  circuit  to  achieve  robust  anti¬ 
tumor  responses  in  our  model  (Figures  7 C,  7D,  and  S6)  and 
that  our  circuit  could  mediate  anti-tumor  effects  when 
delivered  in  vivo  by  multiple  lentiviruses  (Figures  7E,  7F,  and 
S7),  we  hypothesize  that  delivery  efficiencies  will  not  need 
to  be  perfect  to  achieve  a  therapeutic  effect.  However,  we 
note  that  delivery  requirements  will  likely  vary  in  other  tumor 
models. 

Furthermore,  it  will  be  useful  to  test  whether  our  approach 
can  be  used  in  combination  with  other  therapeutic  strategies. 
For  example,  the  efficacy  of  engineered  T  cells  may  be 
enhanced  with  tumor-specific  secretion  of  immunomodulators 
optimized  to  attract  these  cells  to  the  tumor  site  and  activate 
them  while  simultaneously  disrupting  the  immunosuppressive 
microenvironment.  Additionally,  recent  studies  have  shown 
that  many  tumors  express  immunogenic  neoantigens  that  can 
be  targeted  by  the  host  immune  system  (Schumacher  and 
Schreiber,  2015).  We  speculate  that  using  our  strategy  to 
trigger  tumor  killing  in  combination  with  neoantigen-based  im¬ 
munotherapies  could  be  a  fruitful  approach  to  investigate. 
For  example,  our  circuits  could  trigger  robust  killing  of  cancer 
cells  and  release  such  neoantigens.  Furthermore,  secreted 
immunomodulators  could  enhance  neoantigen  presentation 
by  antigen-presenting  cells  and  further  improve  anti-tumor  im¬ 
mune  responses. 

Finally,  since  our  approach  enables  tumor-specific  secretion 
of  defined  combinations  of  immunomodulators,  it  could  be  a 
useful  framework  to  study  the  immunological  mechanisms  that 
underlie  tumor  biology.  Moreover,  synthetic  gene  circuits  have 
been  shown  to  have  the  potential  to  treat  immunological  disor¬ 
ders  by  sensing  extracellular  disease  markers  (Schukur  et  al., 
2015).  Similarly,  our  approach  for  sensing  internal  cell  states 
with  multiple  artificial  promoters  could  be  adapted  to  treat  other 


complex  diseases  that  require  highly  specific  and  multifactorial 
programming  of  immunological  functions  with  simple  modifica¬ 
tions  to  the  synthetic  promoter  inputs  and  immunomodulatory 
outputs. 

In  summary,  modulating  the  immune  system  is  a  promising 
approach  to  treating  complex  diseases.  However,  multiple  dis¬ 
ease-mediated  mechanisms  can  limit  the  effectiveness  of  mono¬ 
therapies  (Mahoney  et  al.,  201 5).  Combinatorial  immunomodula¬ 
tors  can  be  used  to  overcome  this  difficulty,  but  given  the  broad 
effects  of  the  immune  system  on  human  physiology,  it  will  be 
important  to  constrain  the  activity  of  immunomodulators  to  spe¬ 
cific  disease  compartments  to  avoid  undesirable  side  effects 
and  maximize  efficacy.  Thus,  we  envision  that  synthetic  gene 
circuit  immunotherapies  that  are  highly  specific,  tunable,  and 
amenable  to  combinatorial  effector  expression  will  provide 
powerful  strategies  for  engineering  the  immune  system  to  under¬ 
stand  and  treat  disease. 

STAR*METHODS 

Detailed  methods  are  provided  in  the  online  version  of  this  paper 
and  include  the  following: 

•  KEY  RESOURCES  TABLE 

•  CONTACT  FOR  REAGENT  AND  RESOURCE  SHARING 

•  EXPERIMENTAL  MODEL  AND  SUBJECT  DETAILS 

O  Source  and  Culture  of  Primary  Human  T  cells 
O  Animal  Model  Details 
O  Cell  Culture  and  Cell  Lines 

•  METHOD  DETAILS 

o  Plasmid  Construction 
O  Transfection  for  AND  Gate  Optimization 
O  Lentivirus  Production  and  Transduction 
O  Flow  Cytometry 
O  T  Cell-Mediated  Cell  Lysis  In  Vitro 
O  IFN-y  Production  by  T  Cells 

O  Multiple-Output  Circuit  Triggers  STE,  CCL21 ,  IL12,  and 
Anti-PDI  Ab  Production  Specifically  by  OV8  Cells 
O  In  Vivo  Experiments  Validating  Circuit  Therapeutic 
Efficacy 

O  Sponge  Optimization  in  Module  2 
O  miRNA  Backbone  Optimization  in  Module  1 


Figure  6.  Synthetic  Circuits  Triggered  Human  Ovarian  Cancer-Specific  Expression  of  STE,  T  Cell-Mediated  Lysis,  and  IFN-y  Secretion  on 
Tumor  Cells 

(A-C)  The  circuit  in  Figure  5A  was  engineered  to  display  STE  on  cell  surfaces  as  the  output.  STE  was  encoded  in  module  3.  State  “[1,1]  cont.”  indicates  cells 
containing  module  1  and  module  2  and  a  negative  control  module  3con  that  expresses  the  non-specific  transcription  factor  rtTA3  as  the  output.  State  [1,1] 
indicates  cells  containing  module  1 ,  module  2,  and  module  3.  State  [1 ,0]  indicates  cells  containing  module  1 ,  a  negative  control  module  2con  in  which  S(E2F1)p 
expresses  ECFP  without  miRI-BS(B),  and  module  3.  Circuit-transduced  tumor  cells  but  not  normal  cells  expressed  high  levels  of  STE  in  state  [1,1],  which 
triggered  robust  T  cell  killing  and  IFN-y  secretion.  Student’s  t  test  was  performed  to  compare  the  [1 ,1]  states  between  normal  cells  and  OV8  cells,  and  the  [1 ,1] 
state  versus  the  [1 ,0]  state  in  OV8  cells. 

(D-G)  Synthetic  circuits  triggered  human  ovarian  cancer-specific  expression  of  combinatorial  immunomodulators.  The  circuit  in  Figure  5A  was  extended  to 
express  multiple  immunomodulatory  outputs  (surface-displayed  STE  and  secreted  CCL21 ,  IL1 2,  and  anti-PDI  Ab)  specifically  in  OV8  cells  but  not  normal  primary 
cells.  State  [1,1]  cont.  indicates  cells  containing  module  1  and  module  2  and  a  negative  control  module  3con  that  expresses  the  non-specific  TF  rtTA3  as  the 
output.  State  [1,1]  indicates  cells  containing  module  1,  module  2,  and  module  3.  Each  immunomodulatory  gene  was  expressed  from  its  own  G8p  promoter 
encoded  on  a  lentivirus,  except  for  the  anti-PDI  Ab,  which  was  split  into  two  lentiviruses  encoding  LC  and  HC  respectively.  Cell  lines  were  co-infected  with  these 
lentiviral  constructs.  Student’s  t  test  was  performed  to  compare  output  levels  in  the  [1,1]  states  between  normal  cells  and  OV8  cells. 

Error  bars  represent  SEM,  n  =  3  biological  replicates  (*p  <  0.05;  **p  <  0.005). 

See  also  Figure  S4. 
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Figure  7.  Combination  Immunotherapies 
Triggered  by  Cancer-Specific  Circuits 
Significantly  Reduced  Tumor  Burden  and 
Increased  Survival  in  an  Intraperitoneally 
Disseminated  Human  Ovarian  Cancer 
Model  in  Mice 

(A)  OV8  cells  were  transduced  with  constructs 
encoding  a  control  output  rtTA3  (Control),  STE- 
only  (STE),  or  combination  immunomodulators 
(SCIP).  All  groups  were  implanted  into  NSG  mice 
and  injected  periodically  with  human  T  cells.  Stu¬ 
dent’s  t  test  was  performed  to  compare  tumor 
burden  between  groups  at  day  37. 

(B)  Kaplan-Meier  survival  curves  of  various  groups. 
A  log-rank  (Mantel-Cox)  test  was  performed  to 
compare  survival  between  groups. 

(C)  Robust  therapeutic  efficacy  was  achieved  even 
when  only  15%  or  30%  of  tumor  cells  were  trans¬ 
duced  with  the  SCIP.  Student’s  t  test  was  per¬ 
formed  to  compare  tumor  burden  between  groups 
at  day  41 . 

(D)  Kaplan-Meier  survival  curves  of  various  groups 
with  different  percentages  of  the  overall  tumor 
cells  transduced  with  the  SCIP  circuit.  A  log-rank 
(Mantel-Cox)  test  was  performed  to  compare  sur¬ 
vival  between  groups. 

(E)  Lentiviral  delivery  of  the  SCIP-expressing  circuit 
significantly  reduced  ovarian  cancer  burden.  NSG 
mice  were  injected  with  OV8  cells  into  the  perito¬ 
neal  space  at  day  0.  Lentiviruses  encoding  circuit 
module  1  and  module  2  and  a  control  output  rtTA3 
(Control),  or  combination  immunomodulators 
(SCIP)  were  then  injected  i.p.  on  day  7.  These  mice 
were  also  injected  i.p.  periodically  (days  9, 16,  and 
23)  with  human  T  cells.  Student’s  t  test  was  per¬ 
formed  to  compare  tumor  burden  between  groups 
at  day  39. 

(F)  Kaplan-Meier  survival  curves  of  various  groups. 
A  log-rank  (Mantel-Cox)  test  was  performed  to 
compare  survival  between  groups. 

Error  bars  represent  SEM,  n  =  5  biological  repli¬ 
cates  for  all  experiments,  except  n  =  4  biological 
replicates  for  groups  in  (C-D). 

(*p  <  0.05;  **p  <  0.005). 

See  also  Figure  S5,  S6,  and  S7. 
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EXPERIMENTAL  MODEL  AND  SUBJECT  DETAILS 
Source  and  Culture  of  Primary  Human  T  cells 

Human  PBMCs  were  obtained  from  leukoreduction  collar  (Brigham  and  Women’s  hospital  Crimson  Core  Laboratory,  Boston,  MA) 
with  gradient  centrifugation.  To  culture  and  expand  human  T  cells,  human  PBMCs  were  stimulated  with  40  ng/mL  of  anti-human 
CD3  antibody  (clone:  OKT3;  BioLegend,  San  Diego,  CA  ;Catalog  #31 7304)  for  3  days  in  RPM1 1 640  (Invitrogen,  Carlsbad,  CA;  Catalog 
#11875119),  supplemented  with  10%  heat-inactivated  FBS,  10  mM  HEPES  (Life  Technologies  Catalog  #15630080),  0.1  mM  non- 
essential  amino  acids,  1  mM  sodium  pyruvate  (Life  Technologies  Catalog  #11360-070),  1%  Pen/Strep,  50  pM  2-Mercaptoethanol 
(Sigma-Aldrich  Catalog  #M3148-25ML),  and  50  lU/mL  recombinant  human  IL-2  (NCI,  Frederick,  MD). 

Animal  Model  Details 

Animal  studies  were  conducted  at  the  Koch  Institute  Animal  Facility  under  a  protocol  approved  by  MIT’s  Committee  on  Animal  Care. 
Female  NSG  mice  (NOD.Cg-Pr/cc/cSGld  H2rgtm1Wjl/ SzJ;  Stock  No:  005557)  were  purchased  from  the  Jackson  Laboratory  (Bar  Harbor, 
ME).  Mice  were  used  in  experiments  at  the  age  of  8-13  weeks  old. 

Cell  Culture  and  Cell  Lines 

HEK293T,  adult  Human  Dermal  Fibroblast  (aHDF),  MDA-MB-453,  MCF-1 0A,  MCF-1 2A,  and  CCD-841  -CoN  cells  were  obtained  from 
the  American  Type  Culture  Collection,  Rockville,  MD  (HEK293T,  Catalog  #CRL-321 6;  aHDF,  Catalog  #PCS-201  -01 2;  MDA-MB-453, 
Catalog  #HTB-1 31 ;  MCF-1 0A,  Catalog  #CRL-1 031 7;  MCF-1 2A,  Catalog  #CRL-1 0782;  CCD-841  -CoN,  Catalog  #CRL-1 790).  Ovarian 
carcinoma  cell  line  8  (OV8)  engineered  to  stably  express  firefly  luciferase  (OV8-Luc)  was  a  gift  from  Sangeeta  N.  Bhatia  (Massachu¬ 
setts  Institute  of  Technology,  Cambridge,  MA).  Primary  human  ovarian  epithelial  cells  (HOV),  primary  human  ovarian  fibroblasts 
(HOF),  and  primary  human  ovarian  microvasculature  endothelial  cells  (HOMEC)  were  obtained  from  ScienCell  Research  Labora¬ 
tories,  Carlsbad,  CA  (HOV,  Catalog  #7310;  HOF,  Catalog  #7330;  HOMEC,  Catalog  #7300).  Inducible  pluripotent  stem  cell  (iPSC) 
was  a  gift  from  R.  Weiss  (Massachusetts  Institute  of  Technology).  HEK293T,  OV8-Luc,  aHDF,  and  MDA-MB-453  cells  were  cultured 
in  DMEM  (Life  Technologies,  Carlsbad,  CA)  supplemented  with  10%  fetal  bovine  serum  (FBS;  VWR,  Radnor,  PA;  Catalog  #95042- 
108),  1%  Non-Essential  Amino  Acids  (MEM/NEAA;  Hyclone;  Catalog  #16777-186),  and  1%  Pen/Strep  (Life  Technologies  Catalog 
#1 51 40-1 22)  at  37°C  with  5%  C02.  HOV  cells  were  cultured  in  Ovarian  Epithelial  Cell  Medium  (ScienCell  Research  Laboratories  Cat¬ 
alog  #7311)  supplemented  with  Pen/Strep  (ScienCell  Research  Laboratories  Catalog  #0503),  Ovarian  Epithelial  Cell  Growth  Supple¬ 
ment  (ScienCell  Research  Laboratories  Catalog  #7352)  and  were  grown  on  a  poly-L-lysine  (Sigma-Aldrich,  St.  Louis,  MO;  Catalog 
#P4707-50ML)  coated-plate.  HOF  cells  were  cultured  in  Fibroblast  Medium  (ScienCell  Research  Laboratories  Catalog  #2301)  sup¬ 
plemented  with  Pen/Strep  (ScienCell  Research  Laboratories  Catalog  #0503),  Fetal  Bovine  Serum  (ScienCell  Research  Laboratories 
Catalog  #0010),  Fibroblast  Growth  Supplement  (ScienCell  Re-search  Laboratories  Catalog  #2352)  and  were  grown  on  a  poly-L- 
lysine  (Sigma-Aldrich  Catalog  #P4707-50ML)  coated-plate.  HOMEC  cells  were  cultured  in  Endothelial  Cell  Medium  (ScienCell 
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Research  Laboratories  Catalog  #1 001)  supplemented  with  Pen/Strep  (ScienCell  Research  Laboratories  Catalog  #0503),  Fetal  Bovine 
Serum  (ScienCell  Research  Laboratories  Catalog  #0025),  Endothelial  Cell  Growth  Supplement  (ScienCell  Research  Laboratories 
Catalog  #1 052)  and  were  grown  on  a  Fibronectin  (Sigma-Aldrich  Catalog  #F1 1 41 -2MG)  coated-plate.  MCF-1 0A  and  MCF-1 2A  cells 
were  cultured  in  MEGM  BulletKit  (Lonza,  Walkersville,  MD;  Catalog  #CC-3151  and  CC-4136).  All  cell  lines  were  banked  directly  after 
being  purchased  from  vendors  and  used  at  low  passage  numbers.  Cell  lines  were  tested  negative  for  mycoplasma  contamination 
with  IMPACT  PCR  by  the  Division  of  Comparative  Medicine  Diagnostic  Lab  at  MIT. 

METHOD  DETAILS 

Plasmid  Construction 

The  various  circuit  modules,  synthetic  promoters,  and  immunomodulatory  constructs  were  built  using  conventional  restriction 
enzyme  cloning  and  Gibson  assembly.  Plasmids  compositions  used  in  each  experiment  are  detailed  in  supplemental  file  “Plasmids 
by  figures. docx.”  All  DNA  sequences  used  in  this  study  are  detailed  in  GenBank  format  in  the  single  text  file  “All  sequences. docx.” 
Both  files  can  be  downloaded  at  Mendeley  database:  https://data.mendeley.com/datasets/8jsp7jzmt3/draft?a=56eab839-e2b7- 
47d1  -9a6f-ef4d6d1 1 4f97 

Additionally,  individual  GenBank  and  SnapGene  sequence  files  for  the  plasmids  in  this  study  can  be  downloaded  from  a  public 
Dropbox  folder  using  the  link  below: 

https://www.dropbox.com/sh/wbq2y3lx87yzppe/AAAvFM7g6b1  K35Ay1zSreoGDa?dl=0 

Transfection  for  AND  Gate  Optimization 

To  optimize  the  AND  gate  circuit  in  Figures  2,  3,  SI ,  S2,  and  S3,  we  tested  different  combinations  of  two  plasmids:  module  1  variants 
expressed  an  auto-inhibitory  mKate2  RNA  transcript  regulated  by  synthetic  miRNA  (miRI)  and  module  2  variants  expressed  sponges 
for  miRI  within  an  ECFP  transcript,  where  ECFP  enables  us  to  assess  sponge  expression.  In  brief,  1 2  jiL  of  FuGENE  HD  transfection 
reagent  (Promega,  Madison,  Wl;  Catalog  #E2311)  mixed  with  100  jiL  of  OptiMEM  medium  (Life  Technologies  Catalog  #31985)  was 
added  to  the  mixture  of  two  plasmids  (1  pig  each).  During  20  minutes  incubation  of  FuGENE  HD/DNA  complexes  at  room  temperature, 
HEK293T  suspension  cells  were  prepared  and  diluted  to  1 .2  x  1 06  cells/mL  in  culture  medium.  0.5  mL  of  diluted  cells  (0.6  x  1 06  cells) 
were  added  to  each  FuGENE  HD/DNA  complex  tube,  mixed  well,  and  after  5  min  incubation  at  room  temperature  transferred  to  a 
designated  well  in  6-well  plate  containing  2  mL  cell  culture  medium,  followed  by  incubation  at  37°C  with  5%  C02.  Transfected  cells 
were  prepared  for  FACS  analysis  at  48-hour  post-transfection. 

Lentivirus  Production  and  Transduction 

Lentiviruses  were  produced  in  HEK293T  cells  using  co-transfection  in  6-well  plate  format.  In  brief,  12  (iL  of  FuGENE  HD  was  mixed 
with  1 00  fiL  of  OptiMEM  medium  and  was  added  to  a  mixture  of  3  plasmids:  0.5  jig  of  pCMV-VSV-G  vector,  0.5  jag  of  lentiviral  pack¬ 
aging  psPAX2  vector,  and  1  jig  of  lentiviral  expression  vector.  During  20  minutes  incubation  of  FuGENE  HD/DNA  complexes  at  room 
temperature,  HEK293T  suspension  cells  were  prepared  and  diluted  to  3.6  x  106  cells/mL  in  cell  culture  medium.  0.5  mL  of  diluted 
cells  (1 .8  x  1 06  cells)  were  added  to  each  FuGENE  HD/DNA  complex  tube,  mixed  well,  and  incubated  for  5  min  at  room  temperature 
before  being  added  to  a  designated  well  in  6-well  plate  containing  1  mL  cell  culture  medium,  followed  by  incubation  at  37°C  with  5% 
C02.  Media  of  transfected  cells  were  replaced  with  2.5  mL  fresh  culture  media  18  hours  post  transfection.  Supernatant  containing 
newly  produced  viruses  was  collected  at  48-hour  post-transfection,  and  filtered  through  a  0.45  jam  syringe  filter  (Pall  Corporation, 
Ann  Arbor,  Ml;  Catalog  #4614). 

For  testing  of  synthetic  promoters  (Figure  4A  and  4B),  1 :9  dilutions  of  filtered  viral  supernatants  were  prepared  and  3  mL  of  each 
individual  virus  was  used  to  infect  2.5  x  1 05  OV8-Luc,  aHDF  or  HOV  cells  in  the  presence  of  8  jig/mL  polybrene  (Sigma)  overnight.  For 
optimizing  and  testing  the  GAD-expressing  gate  (Figure  5A-5E),  we  transduced  three  viruses  encoding  module  1 ,  module  2,  and 
module  3.  Each  virus  was  diluted  1:3  and  we  used  3  mL  of  the  pooled  viruses  to  infect  2.5  x  105  target  cells  in  the  presence  of 
8  jig/mL  polybrene  (Sigma)  overnight.  Cell  culture  medium  was  replaced  the  next  day  after  infection  and  cells  were  cultured  for  at  least 
one  week  prior  to  FACS  analysis  or  cytotoxicity  assay.  For  creating  cell  lines  bearing  circuits  expressing  multiple  outputs  (Figures  6D- 
6G  and  7),  two  rounds  of  lentiviral  infection  were  performed.  The  viruses  for  the  first  round  consisted  of  module  1 ,  module  2,  and  an 
output  construct  expressing  STE.  The  viruses  for  the  second  round  consisted  of  output  constructs  expressing  all  other  outputs. 

Flow  Cytometry 

For  characterizing  fluorescent  protein  expression,  cells  were  resuspended  with  DM  EM  and  analyzed  by  a  LSRII  Fortessa  cytometer 
(BD  Biosciences,  San  Jose,  CA).  To  characterize  STE  display  and  expression  level  on  cell  surfaces,  various  cell  lines  were  stained 
with  phycoerythrin  labeled  anti-HA  tag  Ab  (Miltenyi  Biotec;  catalog  #130-092-257),  or  anti-His  tag  Ab  (Miltenyi  Biotec;  catalog  #130- 
092-691).  Data  analysis  was  performed  by  FlowJo  software  (TreeStar  Inc,  Ashland,  OR). 

T  Cell-Mediated  Cell  Lysis  In  Vitro 

T  cell-mediated  cell  lysis  was  measured  by  LDH  release  assay  with  CytoTox  96®  Non-Radioactive  Cytotoxicity  Assay  Kit  (Promega, 
Madison,  Wl;  catalog  #G1 780).  Human  T  cells  were  incubated  with  circuit-expressing  target  cells  at  a  effector  to  target  ratio  (E:T  ratio) 
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of  20:1  (6x1  05:3x1  04)  in  duplicate  wells  in  96  well  flat-bottom  plates.  T  cells  and  target  cells  were  also  plated  alone  to  determine  spon¬ 
taneous  LDH  release  from  T  cells  (“effector  spontaneous”)  and  from  target  cells  (“target  spontaneous”).  A  set  of  target  cells  were  also 
plated  and  lysed  with  lysis  buffer  for  determining  maximal  LDH  release  (“target  maximum”).  8.5  hours  later,  cell-free  culture  super¬ 
natant  was  harvested  and  LDH  release  was  determined  following  the  manufacturer’s  protocol.  T  cell-mediated  lysis  was  calculated 
as:  percentage  of  T  cell-mediated  lysis  =  [(experimental  O.D.  value  -  effector  spontaneous  O.D.  value  -  target  spontaneous  O.D. 
value)  /  (target  maximum  O.D.  value  -  target  spontaneous  O.D.  value)]  x  100.  Media  O.D.  values  were  subtracted  from  all  samples 
before  calculation.  STE  mediated  lysis  was  calculated  as:  T  cell  mediated  lysis  of  STE  expressing  cells  -  T  cell  mediated  lysis  of 
non-STE  expressing  cells.  The  results  of  each  experiment  represent  three  biological  replicates. 

IFN-y  Production  by  T  Cells 

Human  T  cells  were  co-cultured  with  tumor  cells  at  an  E:T  ratio  =  20:1  (6x105:3x104)  in  duplicate  wells  in  96  well  flat-bottom  plates. 
Cell-free  medium  was  collected  after  24h  and  IFN-y  concentration  determined  by  Human  IFN-y  DuoSet  ELISA  (R&D  systems,  Min¬ 
neapolis,  MN;  catalog  #DY285). 

Multiple-Output  Circuit  Triggers  STE,  CCL21,  IL12,  and  Anti-PDI  Ab  Production  Specifically  by  OV8  Cells 

To  validate  STE,  CCL21 ,  IL1 2,  and  anti-PDI  Ab  production  by  OV8  cells  but  not  by  normal  cells  (Figure  6D-6G),  2x1 05  OV8-Luc,  HOV, 
or  aHDF  cells  transduced  with  the  SCIP-expressing  circuit  were  plated  in  6-well  plates  and  cell-free  medium  was  collected  after  72h. 
Displayed  STE  on  cell  surfaces  was  measured  via  staining  with  phycoerythrin  labeled  anti-HA  tag  Ab  (Miltenyi  Biotec)  or  phycoery- 
thrin  labeled  anti-His  tag  Ab  (Miltenyi  Biotec)  and  quantified  by  flow  cytometry  analysis.  CCL21  concentration  was  determined  by 
Human  CCL21/6Ckine  DuoSet  ELISA  (R&D  systems;  catalog  #DY366).  IL12  concentration  was  determined  by  Human  IL12  p70 
DuoSet  ELISA  (R&D  systems;  catalog  #DY1270),  and  anti-PDI  Ab  production  was  validated  by  a  competitive  staining  approach. 
The  competitive  staining  was  done  by  incubating  cell-free  medium  (circuit-conditioned  or  control  medium)  with  activated  human 
T  cells  for  30  minutes,  then  followed  by  staining  with  Brilliant  Violet  421  anti-human  CD279  (PD-1)  Ab  [clone:  EH12.2H7;  BioLegend; 
catalog  #329920].  The  percentage  signal  reduction  was  calculated  by  the  following  equation:  (1  -  PD1  staining  MFI  of  T  cells  incu¬ 
bated  with  circuit-conditioned  medium  /  PD1  staining  MFI  of  T  cells  incubated  with  control  medium)  *  100%. 

In  Vivo  Experiments  Validating  Circuit  Therapeutic  Efficacy 

For  in  vivo  experiments,  animals  were  randomly  assigned  to  each  experimental  group.  No  specific  inclusion/exclusion  criteria  were 
used  for  animal  studies.  Investigators  were  not  blinded  to  animal  groups.  Four  to  five  mice  were  used  in  each  experimental  group.  To 
validate  STE  therapeutic  efficacy  on  ovarian  cancer,  OV8-Luc  cells  expressing  doxycycline  (Dox)-inducible  STE  (5x1 05  cells)  were 
injected  intraperitoneally  (i.p.)  into  NSG  mice  on  day  0,  and  STE  expression  were  induced  by  Dox  [0.2mg/mL  Dox  +  2%  Sucrose 
(Sigma-Aldrich;  catalog  #S3929-1KG)  in  drinking  water]  at  day  14.  Activated  human  T  cells  (107  cells)  were  injected  at  day  16 
and  31.  For  measuring  tumor  burden,  mice  were  injected  with  3mg  of  VivoGlo  Luciferin,  In  Vivo  Grade  (Promega;  catalog  #P1043) 
10  minutes  prior  to  bioluminescence  imaging  (IVIS  spectrum  optical  imaging  system,  Xenogen,  Alameda,  CA). 

To  measure  the  efficacy  of  different  therapeutic  output  combinations  on  ovarian  cancer,  OV8-Luc  cells  (5x1 05  cells)  engineered  to 
express  various  outputs  were  injected  i.p.  into  NSG  mice  on  day  0.  Activated  human  T  cells  (107  cells)  were  injected  at  day  7, 14,  21 
and  28.  Tumor  burden  was  monitored  by  IVIS  imaging  using  the  protocol  described  above.  To  measure  survival,  mice  were  eutha¬ 
nized  when  they  developed  signs  of  distress,  such  as  ruffled  fur,  poor  body  posture,  distended  abdomen,  and  jaundice. 

To  identify  the  minimal  percentage  of  tumor  cells  needed  to  be  transduced  with  the  SCIP-expressing  circuit  to  mediate  robust  ther¬ 
apeutic  efficacy,  OV8-Luc  cells  engineered  with  the  SCIP-expressing  circuit  were  mixed  with  wild-type  OV8-Luc  cells  at  various 
ratios  (15%  to  85%;  30%  to  70%).  5x1 05  mixed  OV8-Luc  cells  were  injected  i.p.  into  NSG  mice  on  day  0.  Activated  human 
T  cells  (107  cells)  were  injected  at  day  9,  16,  and  23.  Tumor  burden  was  monitored  by  IVIS  imaging  using  the  protocol  described 
above.  We  used  the  criteria  described  above  to  measure  survival. 

To  deliver  the  SCIP  circuit  with  lentivirus  in  vivo,  lentiviruses  for  each  circuit  component  (module  1 ,  module  2,  STE,  CCL21 +IL1 2, 
and  anti-PDI  Ab,  see  Figure  7E  for  more  details)  were  produced  as  described  above.  For  the  SCI  P  therapy  group,  1 0  mL  of  viruses  for 
each  of  the  five  circuit  components  were  collected  and  concentrated  15-fold  using  Amicon  100K  MW  cutoff  ultrafiltration  columns 
(Millipore,  catalog  #UFC91 0024).  For  the  control  group,  1 0  mL  of  viruses  for  module  1 ,  1 0  mL  of  viruses  for  module  2,  and  30  mL  of 
viruses  encoding  the  control  output  (rtTA3)  was  collected  and  concentrated.  After  concentration,  viruses  were  pooled  and  a  p24 
ELISA  kit  (Advanced  Bioscience  Laboratories,  catalog  #5421)  was  used  to  determine  the  viral  titer.  The  pooled  viruses  had  titers 
of  ~7x  1 010  viral  particles/mL.  Polybrene  was  added  to  the  viruses  to  reach  a  final  concentration  of  1 6  (ig/mL.  600  \.iL  of  viral  super¬ 
natant  was  injected  into  mice  intraperitoneally. 

Sponge  Optimization  in  Module  2 

We  first  optimized  the  miRI  binding-sites  (miRI  -BSs)  in  module  2  for  maximal  sponging  efficiency.  We  designed  1 4  variants  of  miRI  - 
BSs  with  different  mismatches,  insertions,  and  deletions  compared  to  perfect-match  miRI -BSs  (Figures  S1B-S1D).  We  measured 
the  sponging  efficiency  of  each  variant  by  comparing  mKate2  levels  generated  by  module  1  in  the  presence  of  the  module  2  sponge 
encoding  5  identical  tandem  miRI  -BSs  in  the  3'  UTR  of  ECFP  versus  no  miRI  -BSs  (hereafter  reported  as  “fold  ON-OFF  induction”). 
As  an  indication  of  miRI  binding  to  the  sponge,  we  also  expected  to  observe  a  decrease  in  ECFP  expression  in  constructs  where  the 
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3'  UTR  of  ECFP  contained  miRI-BSs  versus  no  miRI-BSs.  We  used  two  highly  active  native  promoters  in  HEK293T  cells,  SSXIp 
(Gure  et  al.,  1997)  and  H2A1p  (Rogakou  et  al.,  1998),  to  express  module  1  and  module  2,  respectively.  In  HEK293T  cells,  H2A1p, 
which  expresses  the  sponge  in  module  2,  is  >  3-fold  stronger  than  SSXIp,  which  expresses  module  1  (Figure  SI  A).  Thus,  the 
module  2  sponge  transcript  concentrations  were  expected  to  be  in  excess  compared  to  the  module  1  transcript,  which  ensured 
enhanced  sponging  efficiency.  When  we  tested  a  sponge  transcript  with  perfect  match  to  miRI  generated  by  module  1 ,  we  found 
a  80%  decrease  in  the  levels  of  ECFP,  even  though  only  a  33%  decrease  was  expected  due  to  the  1 :3  ratio  of  miRI  :sponge  transcript 
concentrations.  This  phenomenon  may  reflect  recycling  of  miRNA  molecules  following  target  transcript  inhibition,  which  is  especially 
efficient  for  perfect-match  binding  sites  (Baccarini  et  al.,  201 1)  (Figure  SI  D). 

We  calculated  the  binding  efficiency  of  miRI  to  its  corresponding  binding  sites  in  the  sponge  (AAG)  with  a  previously  described 
model  (Kertesz  et  al.,  2007),  but  we  found  that  AAG  did  not  predict  sponging  efficiency  (Figure  SIB).  Specifically,  calculations  for 
microRNA  binding  efficiencies  to  target  sites  were  performed  at:  https://genie.weizmann.ac.il/pubs/mir07/mir07_prediction.html. 
For  example,  binding  sites  with  perfect  complementary  to  miRI  (Pe)  had  the  strongest  binding  efficiency  (AAG  =  -26),  but  limited 
mKate2  activation  (Figure  SIC,  1.7-fold  ON-OFF  induction),  which  could  be  a  result  of  efficient  sponge  transcript  degradation 
due  to  perfect  miRNA  complementarity  (Ebert  and  Sharp,  201 0).  A  AAG  =  -15  was  calculated  for  several  binding  sites  that  contained 
various  mismatches  to  miRI  that  resulted  in  bulged  complementarity  to  miRI  (Bv2,  Bv4,  Bv6  and  Bvl  5).  Yet,  the  measured  ON-OFF 
induction  for  these  binding  sites  ranged  from  1 .25-fold  (Bvl  5)  to  2.62-fold  (Bv2),  which  covered  almost  the  entire  dynamic  range  of 
mKate2  activation.  miRI  binding  site  version  Bv2  contained  mismatches  to  nucleotides  9-1 2  of  miRI  and  exhibited  the  best  sponging 
efficiency  (2.6-fold  mKate2  ON-OFF  induction),  likely  because  mismatches  at  these  positions  protect  the  sponge  from  Ago2-medi- 
ated  degradation  while  allowing  efficient  miRI  binding  by  preserving  seed  sequence  complementarity  (Ebert  and  Sharp,  2010).  We 
therefore  used  Bv2  for  the  sponge  in  further  optimization  of  the  circuit.  Hereafter,  we  use  the  abbreviation  ‘#B’  to  refer  to  the  number 
of  tandem  Bv2  binding  sites  separated  by  5  bp  spacers  (for  example,  ‘5B’  refers  to  5  Bv2  binding  sites).  Finally,  we  also  observed  that 
ECFP  levels  could  not  predict  mKate2  activation.  These  results  highlight  the  complexity  of  miRNA-mediated  regulation,  in  which  the 
complementarity  of  the  miRNA  to  its  cognate  binding  sites  determines  not  only  the  inhibitory  pathway  to  which  the  target  mRNA  is 
designated  (mRNA  cleavage,  destabilization  or  translation  inhibition),  but  also  the  rate  of  unloading  of  miRNAs  from  Ago2,  which  is 
enhanced  by  high  complementarity  of  the  miRNA  to  its  target  (De  et  al.,  2013).  While  bulged  BSs  generally  have  better  efficiency  in 
shunting  miRNAs  compared  to  perfect  match  BSs,  we  believe  that  deriving  optimized  rules  on  how  to  introduce  effective  bulges  will 
require  additional  studies. 

We  then  tested  how  the  number,  spacing,  location  and,  architecture  of  miRI  Bv2  binding  sites  in  the  sponge  transcript,  together 
with  sponge  expression  levels,  affect  sponging  efficiency.  We  observed  2.9  and  3.6-fold  mKate2  ON-OFF  induction  when  5B  versus 
10B  binding  sites  were  encoded  in  the  3'  UTR  of  the  module  2  sponge  transcript  expressed  by  H2A1p,  respectively  (Figure  S2A, 
Sponge  1:  H2A1p-ECFP-5B  versus  Sponge  4:  H2A1p-ECFP-10B).  Switching  the  promoter  in  the  sponge  encoding  the  5B  binding 
sites  from  H2A1  p  to  CMVp,  a  4.3-fold  stronger  promoter  compared  to  H2A1  p  (Figure  SI  A),  resulted  in  a  6-fold  ON-OFF  induction  in 
mKate2  levels  (Figure  S2A,  Sponge  3:  CMVp-ECFP-5B).  Therefore,  increasing  the  sponge  transcript  concentration  efficiently  in¬ 
creases  miRI  shunting,  suggesting  that  the  promoter  selected  to  express  module  2  should  be  stronger  than  the  one  selected  to  ex¬ 
press  module  1 .  Switching  the  location  of  the  1 0B  from  the  3'  UTR  to  the  5'  UTR  of  ECFP  in  module  2  (Sponge  4:  FI2A1  p-ECFP-1 0B 
versus  Sponge  5:  FI2A1p-10B-ECFP)  reduced  the  sponging  efficacy  from  3.6-fold  mKate2  ON-OFF  induction  to  2.4-fold  (Figure  S2A), 
which  suggests  that  the  miRI  -BS(B)  were  more  accessible  when  encoded  in  the  3'  UTR  of  the  mRNA  rather  than  in  the  5'  UTR.  Inter¬ 
estingly,  increasing  the  spacing  between  miRI-BSs  located  in  the  3'  UTR  of  the  sponge  transcript  from  5  bp  to  100  bp  (Sponge  2: 
FI2A1  p-ECFP-5B/1  OObp)  only  mildly  reduced  the  sponging  efficiency  from  2.9  to  2.3-fold  mKate2  ON-OFF  induction.  Additionally, 
encoding  24  miRI-BSs  in  the  loop  of  an  RNA  hairpin  with  previously  described  tough-decoy  (TuD)  architecture  (2  miRI-BSs  per 
loop,  12  loops  overall,  Sponge  9:  H2A1  p-ECFP-1 2TuD  (Ebert  and  Sharp,  2010;  Haraguchi  et  al.,  2009)  provided  only  3.8-fold  mKate2 
ON-OFF  induction  (Figure  S2A). 

Finally,  to  increase  the  molar  ratio  of  sponges  produced  from  a  given  promoter,  we  designed  a  module  2  variant  where  10B  were 
encoded  as  an  intron  within  an  ECFP  coding  sequence  in  addition  to  5B  encoded  in  the  3'UTR  of  the  transcript  (Figure  S2A).  In  this 
architecture,  the  donor  and  acceptor  consensus  sequences  that  flank  the  10  intronic  miRI-BS(B)  were  derived  from  a  previously 
described  synthetic  intron  (Sponge  6:  FI2A1  p-ECI -1 0BCons-EC2-5B)  (Greber  et  al.,  2008).  Based  on  ECFP  expression,  this  sponge 
was  correctly  spliced,  thus  producing  ECFP  mRNA  that  contained  5B  and  an  intron  with  10B.  This  sponge  demonstrated  ~4-fold 
mKate2  ON-OFF  induction,  slightly  better  than  Sponge  4:  H2A1  p-ECFP-1 0B. 

We  designed  an  additional  architecture  in  which  the  intronic  consensus  sequences  that  flank  the  10B  were  derived  from  the 
ZRANB1  gene  (Memczak  et  al.,  2013)  (Sponge  7:  H2A1  p-ECI -1  0Bzranb1-EC2-5B).  This  design  was  supposed  to  generate  a  stable 
circular  10B  exon  and  EC1/EC2-5B  introns,  thus  disrupting  ECFP  expression.  This  architecture  demonstrated  enhanced  sponging 
efficacy  with  ~6-fold  mKate2  activation  (Figures  S2A  and  S2B)  and  did  not  produce  ECFP  signal.  To  test  whether  that  the  lack  of 
ECFP  signal  was  caused  by  the  expected  improper  splicing,  resulting  in  an  impaired  ECFP  transcript,  or  because  the  transcript 
was  not  spliced  at  all,  we  generated  an  additional  sponge  in  which  the  ZRANB1  splicing  sequences  were  deleted  (Sponge  8: 
H2A1  p-ECI -1 0B-EC2-5B).  Sponge  8  demonstrated  almost  identical  sponging  efficacy  as  Sponge  7  (Figure  S2A),  which  supports 
the  hypothesis  that  Sponge  7  did  not  splice.  Since  Sponge  7  and  Sponge  8,  which  did  not  splice,  demonstrated  higher  sponging 
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efficiency  than  Sponge  6,  which  did  splice,  we  concluded  that  intronically  encoded  miRNA  sponges  did  not  significantly  contribute  to 
sponging  activity.  Thus,  Sponge  7:  H2A1p-EC1-10BZRANB1-EC2-5B  architecture  was  chosen  as  the  final  ‘optimized  sponge’  design 
for  our  circuit. 

miRNA  Backbone  Optimization  in  Module  1 

We  hypothesized  that  our  circuit  performance  could  be  optimized  by  enhancing  miRI  production  from  module  1 .  We  therefore  modi¬ 
fied  the  pri-miRNA  stems  for  efficient  pri-miRNA  processing  (Figure  S3A  and  Table  S2)  (Auyeung  et  al.,  2013;  Fellmann  et  al.,  2013; 
Suzuki  et  al.,  2015).  Previous  studies  have  described  multiple  features  of  pri-miRNA  stems  that  promote  Drosha-mediated  pri-miRNA 
processing,  including  a  UG  motif  at  the  base  of  the  pri-miRNA  hairpin,  a  narrow  range  of  tolerable  pri-miRNA  stem  lengths,  a  mis¬ 
matched  motif  in  the  basal  stem  region,  and  a  CNNC  motif  16-18bp  downstream  to  the  Drosha  processing  site  that  is  required 
for  efficient  SRp20-dependent  processing  of  the  pri-miRNAs  (Auyeung  et  al.,  2013;  Fellmann  et  al.,  2013).  Our  miRI  design  Mvl 
was  based  on  a  previously  described  intronically  encoded  synthetic  miRNA  (Greber  et  al.,  2008;  Nissim  et  al.,  2014;  Xie  et  al., 
2011)  derived  from  a  miR-30  backbone  (Fellmann  et  al.,  2013).  The  EcoRI  restriction  site  at  position  231-236  was  generated  by 
mutating  the  original  sequence  of  endogenous  miR-30  from  GACTTC  to  GAATTC  (where  position  1  is  the  first  intronic  nucleotide 
following  the  splicing  of  module  1).  This  point  mutation  disrupted  an  essential  CNNC  motif  in  the  miRNA  backbone  and  resulted 
in  reduced  pri-miRNA  processing  efficiency.  The  Xhol  restriction  site  at  position  128-133  was  generated  by  mutating  the  original 
sequence  of  the  endogenous  miR-30  from  CTAAAG  to  CTCGAG,  which  could  further  impair  the  pri-miRNA  secondary  structure. 

We  generated  our  Mv2  design  by  adding  a  CNNC  motif  (CTTCAAGGGGCTA)  downstream  to  base  pair  231  of  Mvl  (Insert  1)  and  an 
additional  insert  (AAGGTATAT)  downstream  to  base  pair  1 33  of  Mvl  (Insert  2)  to  improve  the  secondary  pri-miRNA  structure  that  was 
altered  by  the  formation  of  the  Xhol  restriction  site  (Fellmann  et  al.,  201 3)  (Figure  S3A  and  Table  S2).  To  generate  our  Mv3  design,  we 
added  an  additional  insert  (GACTTC)  downstream  to  base  pair  1 1 1  of  Mvl  (Insert  3)  to  improve  the  secondary  pri-miRNA  structure 
that  was  altered  by  the  formation  of  the  EcoRI  restriction  site  (Fellmann  et  al.,  201 3)  (Figure  S3A  and  Table  S2).  This  resulted  in  a  pri- 
miRNA  backbone  similar  to  the  optimized  miR-E  backbone,  which  is  superior  to  the  miR-30  backbone  (Fellmann  et  al.,  2013). 

We  compared  miRI  production  from  each  backbone  by  characterizing  their  efficiency  at  inhibiting  mKate2  levels  in  the  module  1 
transcript,  in  which  they  were  encoded,  and  ECFP  levels  in  a  module  2  mRNA  containing  5  perfect  match  miRI-BSs  (5Pe)  (Figures 
S3B-S3E).  As  reference,  we  used  a  control  construct  with  no  miRI-BSs  (Figure  S3B).  Pri-miRNA  Mv3  mediated  a  93%  decrease  in 
ECFP  expression  in  the  module  2  construct  with  the  5Pe  sponge  compared  to  the  control  construct,  which  was  superior  to  Mvl  (80% 
inhibition)  and  Mv2  (71%  inhibition),  demonstrating  that  pri-miRNA  design  can  tune  the  efficiency  of  miRI  production.  These  data 
suggest  that  that  miRI  production  was  more  potent  in  Mv3  compared  to  Mv2  and  Mvl . 

We  also  validated  the  effects  of  the  miRI  backbone  on  the  efficiency  of  mKate2  production  by  measuring  mKate2  activation  levels 
in  the  presence  of  the  optimized  sponge  in  module  2  (Figures  S3D  and  S3E).  The  absolute  mKate2  signal  in  the  OFF  and  ON  states 
(e.g.,  mKate2  levels  measured  with  control  versus  optimal  sponge)  were  both  affected  by  the  pri-miRNA  design:  Mvl  exhibited  an 
OFF-ON  shift  from  ~1 400  to  ~8400  Mean  Fluorescence  Intensity  (MFI)  (5.8-fold  activation  ratio),  Mv2  from  ~1 70  to  ~500  MFI  (3-fold 
activation)  and  Mv3  from  ~380  to  ~2600  MFI  (6.8-fold  activation).  Based  on  the  ECFP  repression  data  from  comparing  the  5Pe 
sponge  to  the  control  construct,  we  expected  Mv3  to  be  more  efficient  at  inhibiting  mKate2  than  Mvl  and  Mv2  in  the  presence  of 
the  5Pe  sponge.  However,  Mv2  exhibited  potent  inhibition  of  mKate2  levels  in  both  the  OFF  and  ON  states.  This  suggests  that 
our  modifications  to  the  pri-miRNA  design  may  not  only  affect  the  efficiency  of  miRI  production,  but  also  other  regulatory  mecha¬ 
nisms.  We  hypothesize  that  these  pri-miRNA  modifications  could  modify  the  efficiency  of  existing  splicing  sites  or  may  have  added 
alternative  splicing  sites  in  the  miRNA  intron  that  could  impair  mKate2  production  by  alternative  splicing.  This  hypothesis  is  supported 
by  the  fact  that  the  pri-miRNAs  in  designs  Mv2  and  Mv3  include  a  potential  splicing  donor  sequence  AAGGT  introduced  by  Insert  2, 
and  that  a  potential  donor  site  CAGGT  in  Mvl  and  Mv2  was  modified  to  CAGGG  by  Insert  3  in  Mv3  (Table  S2). 

As  noted  in  the  main  text,  we  selected  the  Mv3  architecture  for  further  experimentation  since  it  had  low  mKate2  level  in  OFF  state 
and  exhibited  the  highest  fold  activation  by  the  sponge  in  ON  state. 

Finally,  the  miRNA  duplex  sequence  can  be  modified  to  promote  asymmetric  guide  strand  selectivity  that  favors  miRI  production. 
The  miRI  pri-miRNA  is  processed  into  a  miRNA  duplex  capable  producing  both  miRI  encoded  by  the  duplex  3p  arm  and  an  alter¬ 
native  miRNA  encoded  by  the  duplex  5p  arm  (Table  S2).  The  ratio  of  3p  (miR1):5p  production  is  determined  by  the  superposed  pat¬ 
terns  of  the  5'  end  nucleotide  identity  and  the  thermodynamic  stability  of  two  miRNA  duplex  termini  (Suzuki  et  al.,  2015).  Quantitative 
thermodynamic  and  sequence  analysis  of  the  miRI  duplex  has  shown  that  the  3p  arm  has  5'-uridine  and  a  thermodynamically 
unstable  5'  end,  whereas  5p  arm  has  5'-guanine  and  a  thermodynamically  stable  5'  end  (Suzuki  et  al.,  2015).  This  indicates  that 
the  miRI  duplex  is  optimized  to  produce  mostly  the  3p  arm  (e.g.,  miRI);  thus,  it  was  therefore  not  modified.  Nevertheless,  we 
note  that  asymmetric  guide  strand  selectivity  is  an  important  consideration  in  designing  synthetic  miRNAs. 

Bioinformatics  Framework  for  Synthetic  Promoters  Design 

We  identified  genes  that  were  differentially  expressed  in  human  ovarian  cancer  cells  compared  to  healthy  tissues  from  the  publicly 
available  The  Cancer  Genome  Atlas  (TOGA)  database  (Bell  et  al.,  2011)  using  R  (version  3.0.2).  We  followed  standard  procedures  to 
calculate  differentially  expressed  genes  between  cancer  cases  compared  to  controls  using  t  test.  The  p  values  were  corrected  for 
multiple  testing  using  the  Bonferroni  correction.  We  then  selected  genes  that  were  previously  identified  as  transcription  factors 
(TFs)  from  the  literature  (Vaquerizas  et  al.,  2009;  Wingender  et  al.,  2013).  From  these  genes,  we  selected  the  ones  that  were 


e6  Cell  171,  1 138-1 150.e1-e7,  November  16,  2017 


Cell 


hypomethylated  in  cancer.  For  that  purpose,  we  used  the  publicly  available  database  of  DNA  Methylation  and  gene  expression  in 
Human  Cancer  MethHC  (Huang  et  al.,  2015).  This  database  provides  the  most  hyper-  and  hypomethylated  genes  for  different  can¬ 
cers.  We  selected  TF  genes  with  the  most  hypomethylated  promoter  regions  using  default  settings.  TF  expression  in  OV8  cells  was 
confirmed  using  the  CCLE  database  (Barretina  et  al.,  2012).  We  then  determined  the  binding  motifs  for  ovarian  cancer-enriched  TFs 
using  GREAT,  MEME,  JASPAR  and  MOTIFMAP  (Bailey  et  al.,  2009;  Daily  et  al.,  201 1 ;  McLean  et  al.,  201 0;  Sandelin  et  al.,  2004)  and 
assembled  them  into  our  synthetic  promoters  (Figures  4A  and  4B). 

Analysis  of  Natural  and  Synthetic  Promoters 

We  analyzed  promoter  sequences  using  Matlnspector  software  to  predict  transcription  factor  binding  sites  (TF-BSs)  (Cartharius 
et  al.,  2005).  We  ran  our  analysis  on  two  natural  promoters  (SSX1  p,  H2A1  p)  and  two  synthetic  promoters  (S(E2F1)p,  S(cMYC)p)  using 
core  similarity  0.9  and  matrix  similarity  optimized  +  0.05.  We  examined  individual  TF-BSs  as  well  as  TF-BS  families  (defined  as  TFs 
with  highly  similar  TF-BS  -  for  example,  the  binding  sites  for  MYC,  USF  and  Max  were  grouped  as  the  family  E-BOX).  We  found  that 
the  synthetic  promoters  had  more  predicted  TF-BSs  per  1 0Obp.  For  example,  S(cMYC)p  contained  35.8  TF-BS  matches  per  1 00  bp 
while  H2A1p  and  SSXIp  contained  ~7.9  TF-BS  matches  per  lOObp.  The  synthetic  promoters  had  matches  to  less  TF-BS  families 
(~6-7  different  TF-BS  families)  compared  with  the  natural  promoters  (~40-41  different  TF-BS  families).  The  synthetic  promoters 
had  more  TF-BS  matches  for  each  of  the  different  TF-BS  families  compared  with  the  natural  promoters.  Furthermore,  we  identified 
the  TF-BS  family  with  the  highest  number  of  matches  for  each  promoter,  and  found  that  the  top  TF-BS  family  for  the  synthetic 
promoters  had  more  TF-BS  matches  than  the  natural  promoters.  For  example,  E2F  TF-BS  family  was  found  to  be  66%  of  the  total 
predicted  TF-BSs  in  S(E2F1)p  and  the  E-BOX  TF-BS  family  (corresponding  to  MYC,  USF,  MAX)  were  found  to  be  29%  of  the  total 
predicted  TF-BSs  in  S(cMYC)p.  In  contrast,  the  ZF5F  (ZFP  1 61 )  TF-BS  family  was  found  to  be  only  1 1  %  of  the  total  predicted  TF-BSs 
in  H2A1  p  while  the  MYBL  (MYB)  TF-BS  family  was  found  to  be  only  6%  of  the  total  predicted  TF-BSs  in  SSX1  p. 

Multi-output  AND  Gate  Optimization  for  Combinatorial  Immunomodulation 

We  first  characterized  the  output  levels  of  four  distinct  circuit  states  in  OV8  cells  (Figure  5B).  To  do  so,  we  encoded  each  module  on 
separate  lentiviruses  and  infected  OV8  cells  with  the  corresponding  module  3  architecture  and  either:  (A)  a  negative  control  module  1 
that  expresses  the  non-specific  transcription  factor  rtTA3  (module  Icon)  and  a  negative  control  module  2  in  which  S(E2F1)p  ex¬ 
presses  ECFP  without  miRI-BS(B)  (module  2con)  to  model  the  scenario  when  neither  cancer-specific  promoters  are  active  [0,0]; 
(B)  module  Icon  and  module  2  to  model  the  scenario  when  only  S(E2F1)p  is  active  [0,1];  (C)  module  1  and  module  2con  to  model 
the  scenario  when  only  S(cMYC)p  is  active  [1,0];  (D)  module  1  and  module  2  to  model  the  scenario  when  both  promoters  are 
active  [1,1]. 

The  circuits  expressed  minimal  mKate2  fluorescence  in  states  [0,0]  and  [0,1].  In  state  [1 ,0],  there  was  a  low  background  mKate2 
level  that  increased  from  5  to  8  GAL4-BSs  in  the  synthetic  GALp  and  was  attenuated  with  the  presence  of  miRI  -BS(Pe)  in  module  3 
(Figure  5B).  Only  when  both  input  constructs  were  present  (state  [1,1])  did  we  observe  a  significantly  increased  mKate2  output. 
Furthermore,  the  mKate2  output  level  could  be  enhanced  by  increasing  the  number  of  GAL4-BSs  in  the  module  3  synthetic  promoter 
from  5  to  8  GAL4-BSs.  Although  adding  perfect  match  miRI  binding  sites  in  the  3'  UTR  of  the  mKate2  transcript  decreased  the  basal 
and  maximal  output  levels,  it  also  improved  fold  activation  (G5-Pe,  G8-Pe,  and  G14-Pe  versus  G5,  G8,  and  G14,  respectively).  The 
G8-Pe  circuit  exhibited  negligible  mKate2  levels  in  states  [0,0]  and  [0,1],  and  high  8-fold  mKate2  activation  in  state  [1,1]  over  the 
minimal  mKate2  expression  in  state  [1 ,0],  and  was  thus  chosen  for  further  use. 

QUANTIFICATION  AND  STATISTICAL  ANALYSIS 

Statistical  analysis  was  done  by  Prism  software  version  7.01  (GraphPad  Software,  Inc.,  La  Jolla,  CA). 

DATA  AND  SOFTWARE  AVAILABILITY 

Plasmids  compositions  used  in  each  experiment  are  detailed  in  supplemental  file  “Plasmids  by  figures. docx”.  All  DNA  sequences 
used  in  this  study  are  detailed  in  GenBank  format  in  the  single  text  file  “All  sequences. docx”.  Both  files  can  be  downloaded  at  Men- 
deley  Data:  https://data.mendeley.com/datasets/8jsp7jzmt3/draft?a=56eab839-e2b7-47d1  -9a6f-ef4d6d1 1 4f97. 

Additionally,  individual  GenBank  and  SnapGene  sequence  files  for  the  plasmids  in  this  study  can  be  downloaded  from  a  public 
Dropbox  folder  using  this  link:  https://www.dropbox.com/sh/wbq2y3lx87yzppe/AAAvFM7g6b1K35Ay1zSreoGDa7dk0. 
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Figure  SI.  Optimization  of  the  miRNA  Binding-Site  Sequences  in  the  Sponge  in  Module  2,  Related  to  Figure  3 

(A)  The  activity  of  the  promoters  SSXIp,  H2A1p,  and  CMVp  in  HEK293T  cells. 

(B-C)  Various  miRNA  binding  sequences  (miRI  -BSs)  used  in  the  sponge  affect  sponging  activity,  assayed  by  mKate2  fold  activation  by  each  sponge  (state  [1,1] 
versus  state  [1 ,0]  in  Figure  2B).  Nucleotides  highlighted  in  dark  blue  are  mismatched  nucleotides.  Nucleotides  highlighted  in  red  are  insertions.  Light  blue 
capital-delta  (A)  indicates  deleted  nucleotides.  SSXIp  was  used  to  drive  module  1  expression  and  H2A1p  was  used  to  drive  module  2  expression.  Module  1 
encoded  a  transcript  with  miRI  as  an  intron  within  the  mKate2  gene,  where  miRI  targeted  3  perfect  match  miRI  -BS  encoded  downstream  of  mKate2.  Module  2 
encoded  the  miRI  sponges,  described  above,  downstream  of  an  ECFP  gene. 

(D)  MFI  of  mKate2  and  ECFP  in  binding  site  optimization  experiments.  Error  bars  represent  SEM,  n  =  3  biological  replicates.  MFI  =  Mean  Fluorescence  Intensity. 
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Figure  S2.  Optimization  of  the  Sponge  Architecture  in  Module  2,  Related  to  Figure  3 

(A)  The  number  of  miRNA  binding  sites,  the  sponge  expression  level,  and  the  architecture  of  the  sponge  affect  sponging  activity.  Fold  activation  was  assayed  by 
mKate2  fold  activation  by  each  sponge  (state  [1 ,1]  versus  state  [1 ,0]  in  Figure  2B).  SSX1  p  was  used  to  drive  module  1  expression  and  H2A1  p  or  CMVp  was  used  to 
drive  module  2  expression,  as  indicated.  Module  1  encoded  a  transcript  with  miRI  as  an  intron  within  the  mKate2  gene,  where  miRI  targeted  3  perfect  match 
miRI  -BS  encoded  downstream  of  mKate2.  Module  2  encoded  the  miRI  sponges.  5B  and  1 0B  refers  to  5  or  1 0  bulged  miRI  binding  sites  in  the  sponge  based  on 
the  Bv2  design  from  Figure  SI  B.  EC1  and  EC2  correspond  to  exon  1  and  exon  2  of  ECFP  and  constructs  in  which  1 0B  was  encoded  within  the  ECFP  transcript 
with  different  architectures  (Cons  and  ZRANB1).  1 2TuD  (Tough  Decoy)  architecture  refers  to  24  miRI  -BSs  encoded  within  1 2  loop  of  an  RNA  hairpin  (2  miRI  -BSs 
per  hairpin  loop).  5B/100  bp  refers  to  5  bulged  miRI -BSs  separated  by  lOObp  linkers  rather  than  5bp  in  all  other  sponges.  Error  bars  represent  SEM,  n  =  3 
biological  replicates. 

(B)  MFI  of  mKate2  and  ECFP  in  sponge  architecture  optimization  experiments. 
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Figure  S3.  Pri-miRNA  Backbones  Used  to  Produce  MiRI  from  Module  1  Affect  Circuit  Performance,  Related  to  Figure  3 

(A)  Module  1  variants  Mv1-Mv3  expressed  a  transcript  in  which  the  different  miRI  pri-miRNA  architectures  shown  were  encoded  as  an  intron  within  the  mKate2 
gene,  where  miRI  then  targeted  3  perfect  match  miRI  -BS  encoded  downstream  of  mKate2.  The  colored  circles  denote  the  location  of  the  insertion  sequences  for 
creating  Mv2  and  Mv3  (See  Table  S2  for  detailed  sequences  and  functional  pri-miRNA  modifications). 

(B-E)  Different  pri-miRNA  backbones  modulate  the  production  of  miRNA  and  the  ON:OFF  ratio  of  the  circuit.  (B-C)  We  tested  the  three  different  module  1  variants 
together  with  module  2  variations  composed  of  a  control  sponge  with  no  miRI  binding  sites  (‘Cont.’)  and  a  sponge  with  5  perfect  match  miRI  binding  sites 
downstream  to  an  ECFP  gene  (‘5Pe’).  Module  1  was  expressed  by  SSX1  p  and  module  2  was  expressed  by  H2A1  p.  Mvl  -Mv3  exhibited  different  miRI  production 
efficiencies,  as  measured  by  ECFP  signal  reduction  with  the  5Pe  module  2  compared  to  ECFP  signal  in  the  control  module  2;  Mv3  had  the  most  efficient  miRI 
production. 

(D-E)  mKate2  expression  levels  with  control  sponge  or  the  optimized  sponge  (‘Sponge  7’,  Sponge  7  from  Figure  S2A).  Mv3  exhibited  a  low  mKate2  level  in  the 
OFF  state  and  the  highest  fold  activation  by  the  sponge  in  the  ON  state. 

Error  bars  represent  SEM,  n  =  3  biological  replicates. 
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Figure  S4.  Multi-output  Circuit  Specifically  Triggered  T  Cells  to  Kill  Ovarian  Tumor  Cells  and  Secrete  IFN-y,  Related  to  Figures  5  and  6 

(A)  FACS  histograms  showing  the  activity  of  single  promoters  and  the  G8-Pe  circuit  when  transduced  into  T  cells  via  lentivirus.  The  percentages  on  the  histograms 
denote  the  percent  of  cells  with  higher  mKate2  expression  level  than  the  value  of  the  dotted  line.  A  constitutive  human  ubiquitin  C  promoter  expressing  mKate2 
was  used  for  the  “infection  control”  group. 

(B)  S(cMYC)p,  S(E2F1)p,  and  the  G8-Pe  circuit  exhibit  high  activity  in  OV8  cells  but  not  in  various  non-tumorigenic  cell  lines.  OV8  is  an  ovarian  cancer  cell  line.  FIOF 
(human  ovarian  fibroblasts),  HOMEC  (human  ovarian  microvasculature  epithelial  cells),  CCD-841 -CoN  (non-tumorigenic  colonic  epithelial  cells),  and  MCF-10A 
and  MCF-12A  (non-tumorigenic  mammary  epithelial  cells),  are  non-tumorigenic  cell  lines.  iPSCs  =  inducible  pluripotent  stem  cells. 


(legend  continued  on  next  page) 
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(C-E)  T  cell-mediated  killing  of  G8-Pe  circuit-transduced  tumor  cells  was  significantly  greater  than  of  normal  cells.  The  percentage  of  cell  death  mediated  by 
T  cells  from  donor  #1  (D),  donor  #2  (E)  and  the  calculated  mean  (C)  are  shown.  *  indicates  condition  not  tested.  Note  that  total  cell  death  and  STE-mediated  cell 
death  of  the  [1 ,0]  state  in  aHDF  and  HOV  cells  were  not  measured  because  this  state  simulates  a  condition  in  which  only  one  of  the  input  promoters  is  active  and  is 
expected  to  result  in  lower  killing  than  state  [1,1],  in  which  both  input  promotes  are  active.  Total  cell  death  and  STE-mediated  cell  death  correspond  to  un¬ 
corrected  cell  death  and  background-killing  corrected  cell  death,  respectively.  See  STAR  Methods  for  more  details. 

(F)  G8-Pe  circuit-transduced  tumor  cells  trigger  robust  IFN-y  secretion  by  T  cells.  Data  shown  represent  IFN-y  levels  secreted  by  T  cells  from  donor  #1  and 
donor  #2.  Error  bars  represent  SEM,  n  =  3  biological  replicates. 

(G-l)  aHDF,  HOV,  and  OV8  constitutively  displaying  STE  can  trigger  robust  STE-mediated  killing  and  IFN-y  secretion  by  T  cells.  A  constitutive  human 
ubiquitin  C  promoter  (hubCp)  was  used  to  drive  STE  expression  on  aHDF,  HOV,  and  OV8.  All  three  cell  lines  displayed  STE  and  triggered  strong  STE-mediated 
cell  killing  and  IFN-y  secretion  by  T  cells.  Error  bars  represent  SEM,  n  =  3  biological  replicates. 

(J)  Doxycycline  (Dox)-inducible  expression  of  STE  triggered  significant  T  cell-mediated  killing  in  vivo.  Inducible  display  of  STE  with  doxycycline  (Dox)  on  the 
surface  of  OV8  cells  triggered  human  T  cell  mediated  killing  in  vivo  in  a  mouse  model  of  disseminated  intraperitoneal  ovarian  cancer.  Anti-tumor  activity  was 
abrogated  when  STEs  were  not  displayed  (T  cells  only)  or  T  cells  were  not  injected  (Dox  only).  Student’s  t  test  was  used  to  compare  tumor  burdens  between 
groups  at  day  45.  Error  bars  represent  SEM,  n  =  5  biological  replicates. 

(K)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  the  Dox  +  T  cell  group  in  Figure  S4J. 

(L)  OV8  cells  engineered  with  a  Dox-inducible  STE  construct  expressed  high  levels  of  STE  on  their  cell  surfaces  after  induction.  Histograms  of  STE  staining  on  OV8 
cells  are  shown.  The  gray  and  orange  histograms  represent  the  staining  of  uninduced  OV8  and  Dox-induced  OV8  cells,  respectively.  Data  shown  are  repre¬ 
sentative  histograms,  which  were  repeated  in  3  independent  experiments. 

(*p  <  0.05;  **p  <  0.005;  ***p  <  0.001). 
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Figure  S5.  Combination  Immunotherapies  Triggered  by  the  Circuits  Significantly  Reduced  Tumor  Burden  in  an  Intraperitoneally  Dissemi¬ 
nated  Ovarian  Cancer  Model  within  NSG  Mice  Periodically  Injected  with  Human  T  Cells,  Related  to  Figure  7 

(A)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  each  tested  group  (n  =  5  biological  replicates).  X-axes  denote  days  post  tumor  injection,  Y-axes 
denote  tumor  burden  (photons/sec). 

(B)  In  vivo  tumor  bioluminescence  images  of  each  mouse  from  day  5  (D5)  to  37  (D37).  These  data  correspond  to  the  experiment  shown  in  Figure  7A  and  7B. 
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Figure  S6.  Robust  Therapeutic  Efficacy  Was  Achieved  Even  When  Only  15%  of  Tumor  Cells  Were  Transduced  with  the  SCIP-Expressing 
Circuit  and  Then  Implanted  into  the  Peritoneal  Cavity  of  Mice,  Related  to  Figure  7 

(A)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  each  tested  group  (n  =  4  biological  replicates).  X-axes  denote  days  post  tumor  injection,  Y-axes 
denote  tumor  burden  (photons/sec). 

(B)  In  vivo  tumor  bioluminescence  images  of  each  mouse  from  day  5  (D5)  to  41  (D41).  These  data  correspond  to  the  experiment  shown  in  Figure  7C  and  7D. 


Cell 


SCIP: 


S(cMYC)p' 

S(E2F1)p. 


^  STE  IPe 

G8p 

_  j^lU2P2^C^11Pe 
G8p 

r+  a 

G8p 


Days  post  tumor  injection  Days  post  tumor  injection 


Figure  S7.  Lentiviral  Delivery  of  SCIP  Circuit  Mediated  Robust  Therapeutic  Efficacy  in  an  Intraperitoneally  Disseminated  Human  Ovarian 
Cancer  Model  in  Mice,  Related  to  Figure  7 

(A)  Tumor  burden  of  individual  mice  (black)  and  mean  (red)  of  each  tested  group  (n  =  5  biological  replicates).  X-axes  denote  days  post  tumor  injection,  Y-axes 
denote  tumor  burden  (photons/sec). 

(B)  In  vivo  tumor  bioluminescence  images  of  each  mouse  from  day  13  (D13)  to  39  (D39).  These  data  correspond  to  the  experiment  shown  in  Figure  7E  and  7F. 


